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ABSTRACT 

We compare line emission calculated from theoretical disk models with optical to 
sub-millimeter wavelength observational data of the gas disk surrounding TW Hya and 
infer the spatial distribution of mass in the gas disk. The model disk that best matches 
observations has a gas mass ranging from ~ 10~^ — 10~^Mq for 0.06AU < r < 3.5AU 
and ~ O.O6M0 for 3.5AU < r < 200AU. We find that the inner dust hole (r < 3.5AU) 
in the disk must be depleted of gas by ~ 1 — 2 orders of magnitude compared to 
the extrapolated surface density distribution of the outer disk. Grain growth alone is 
therefore not a viable explanation for the dust hole. CO vibrational emission arises 
within r ~ 0.5AU from thermal excitation of gas. [01] 6300A and 5577A forbidden 
lines and OH mid-infrared emission are mainly due to prompt emission following UV 
photodissociation of OH and water at r < 0.1 AU and at r ~ 4AU. [Nell] emission 
is consistent with an origin in X-ray heated neutral gas at r < 10 AU, and may not 
require the presence of a significant EUV {hu > 13.6eV) flux from TW Hya. H2 pure 
rotational line emission comes primarily from r ~ 1 — 30AU. [OI]63/im, HCO"*" and CO 
pure rotational lines all arise from the outer disk at r ~ 30 — 120AU. We discuss planet 
formation and photoevaporation as causes for the decrease in surface density of gas and 
dust inside 4 AU. If a planet is present, our results suggest a planet mass ~ 4 — 7Mj 
situated at ~ 3AU. Using our photoevaporation models and the best surface density 
profile match to observations, we estimate a current photoevaporative mass loss rate of 
4 X 10~^Mq yr~^ and a remaining disk lifetime of ~ 5 million years. 
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Subject headings: Astrochemistry — line: formation — protoplanetary disks — planet- 
disk interactions — stars:individual (TW Hya) 

1. Introduction 

Circumstellar disks are widely believed to go through a transition epoch as they evolve from 
an initial, gas-rich, classical T Tauri (CTTS) phase to the weak-lined T Tauri (WTTS) debris disk 
phase (e.g., Strom et al. 1989, Skrutskie et al. 1990, Najita et al. 2007). Transition disks, which 
typically are a small fraction (~10-20%) of all disks, are characterized by an outer, optically thick 
dust disk and an inner dust-depleted region or "hole". However, it is unclear if all disks go through 
the transition phase (e.g., Currie et al. 2009, Muzerolle et al. 2010). Transition disks are sometimes 
found to be accreting at rates similar to young disks suggesting that the dust holes are not devoid 
of gas, but many objects do not show signs of accretion (Sicilia-Aguilar, Henning & Hartmann 
2010). The duration of the transition epoch is also uncertain (e.g., Muzerolle et al. 2010), although 
a recent study of the statistics of spectral energy distributions suggests times of ~ a few 10^ years 
(Ercolano &: Clarke 2010). Studying these partially evacuated disks can potentially provide many 
insights into how disks evolve and form planets. 

At least three different scenarios have been proposed to explain the morphology of transition 
disks (see Najita et al. 2007 for a summary). First, the infrared opacity hole may not be a true 
depletion of material, but the result of dust evolution. Coagulation and grain growth in disks might 
lead to the formation of larger, "invisible" solid objects, thereby emptying the inner regions of small, 
infrared-emitting dust. In this sequence of events, gas and large solids or planetesimals remain and 
there is no true mass depletion (e.g., Dullemond & Dominik 2005). In a second situation, the hole 
in the disk is caused by dynamical interactions with an embedded planet (e.g., Skrutskie et al. 
1990, Marsh & Mahoney 1992, Takeuchi et al. 1996, Varniere et al. 2006). The Jovian or larger 
mass planet first creates a gap in the disk and then accretes most of the accretion flow from the 
outer disk, while a small fraction continues past the planet to accrete onto the star (e.g., Lubow 
& D'Angelo 2006). The third explanation for the inner hole is photoevaporation, a process by 
which high energy radiation from the central star heats the gas to escape speeds and creates gaps 
(and then quickly central holes) a few AU from the star (e.g., Alexander et al. 2006a, 2006b, 
Gorti, Dullemond & Hollenbach 2009, Owen et al. 2010). Determining which of the above three 
mechanisms dominate in the creation of the inner opacity hole may be key to deciphering the role 
played by transition disks in disk evolution, planet formation, and disk dispersal. 

TW Hya is the closest known (at 51 pc, Mamajek 2005) classical T Tauri star with a nearly 
face-on (e.g., Qi et al. 2006) transition disk that exhibits an inner hole in its (IR-emitting) dust 
distribution (Calvet et al. 2002, Ratzka et al. 2007, Hughes et al. 2009). It is a very well-studied 
object - the central star and disk have been investigated from the X-ray to centimeter wavelengths. 
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The X-ray, Far Ultraviolet (FUV; 6eV< hu <13.6eV) and optical flux from the central star are 
well known (e.g., Kastner et al. 2002, Herczeg et al. 2002, Stelzer &; Schmitt 2004, Robrade & 
Schmidt 2006). The duslQ disk has been studied and extensively modeled, with one of the earliest 
detailed models being that of Calvet et al. (2002, hereafter C02) who inferred the presence of a dust 
depleted inner cavity of size ~ 4AU from the spectral energy distribution. A very small amount of 
dust persists in the inner hole, and the observed Wfiui silicate emission feature has been attributed 
to this component (C02). Ratzka et al. (2007) also model the dust disk to infer a dust cavity, 
but with a smaller size, ~ 0.5 — 0.8 AU. The presence of a dust hole, or at least a significantly 
dust depleted inner cavity, has been confirmed by imaging and interferometric studies (Eisner et 
al. 2006, Hughes et al. 2007). The star has a low, but measurable, accretion rate with estimates 
ranging from 4 x 10~^° Mq yr"^ (Muzerolle et al. 2000) to > 10"^ Mq yr"^ (Alencar & Batalha 
2002, Eisner et al. 2010), indicating the presence of significant gas in the inner dust hole. 

The nature of the transition disk around TW Hya is uncertain. While a dust hole is clearly 
present and gas must be present in the hole to account for the observed accretion, the amount of 
gas in this region is not well determined. A radial velocity planet was reported by Setiawan et al. 
(2008), but located closer to the star (~ 0.04AU) than a planet that could cause the 4 AU dust 
edge. (This detection has been disputed and the radial velocity variations attributed to stellar 
surface activity, see Huelamo et al. 2008, Figueira et al. 2010.) TW Hya is known to be a strong 
far ultraviolet (FUV) and X-ray source and observed to have outward flows traced by [NeII]12.8^m 
emission (Herczeg et al. 2007, Pascucci & Sterzik 2009), which implies that the disk is presently 
undergoing photoevaporation. However, photoevaporation theory has traditionally predicted a 
clean evacuated inner region devoid of dust and gas soon after the formation of a hole, which the 
observed accretion rate onto the central star contraindicates. Owen et al. (2011) argue that the 
inner disk draining epoch, although short (< 10^ years, the local viscous timescale at 1 — 4 AU 
where the gap opens), may still be observed as an accreting transition disk. 

Clearly, spectral line observations of gas in the TW Hya disk might distinguish between the 
above scenarios. In recent years, many detections of gas emission lines from TW Hya have become 
available. The disk around TW Hya was detected in the pure rotational transitions of CO by 
Zuckerman et al. (1995), van Zadelhoff et al. (2001) and later by Qi et al. (2004, 2006) who 
also imaged the disk. Herczeg et al. (2002) detected UV emission from H2 from the disk surface 
near the star. Emission in the near infrared from vibrationally excited H2 was detected by Bary 
et al. (2003). Thi et al. (2004) detected sub-millimeter emission from the molecules CO, HCO+, 
DCO+, CN and HCN. CO rovibrational emission from the disk has been observed (Rettig et al. 
2004, Salyk et al. 2007, Pontoppidan et al. 2008). Resolved line emission of [Nen]12.8^m was 
reported by Herczeg et al. (2007) and Pascucci & Sterzik (2009). Bitner et al. (2008), as part 
of a ground-based survey using TEXES (Richter et al. 2003), failed to detect pure rotational H2 



^We refer to solid particles with sizes smaller <lmm as "dust". Larger solids are ignored as they have too little 
surface area and are therefore not important for gas physics and topics relevant to this paper. 



- 4 - 



emission from the disk, although Najita et al. (2010) detected emission with the larger beam of 
Spitzer. Najita et al. (2010) present a high-quality Spitzer IRS spectrum with numerous emission 
features of Ha, OH, CO2, HCO+, [Nell], [Nelll] and HI. Thi et al. (2010) recently reported the 
Herschel PACS detection of the [OI]63/im line from TW Hya. In addition to this rich dataset, 
strong optical emission has been detected from [01] 6300 A and 5577 A (Al encar &: Batalha 2002, 
S. Edwards, priv. comm.). Pascucci et al. (2011) reanalyze the optical data to also report the 
detection of emission from [SII] 4069 A and [Mgl 4571 A. 

The abundant observational data available for the TW Hya system make it an excellent target 
for studying the physical nature of transition disks through theoretical modeling. Gas emission 
line strengths are very sensitive to the density, temperature and the abundance of the emitting 
species (e.g., Gorti & Hollenbach 2008). The emission lines observed largely originate from trace 
species, and in order to interpret the data accurately, detailed models of gas chemistry and emission 
are needed. In this paper, we use our previous disk models (Gorti & Hollenbach 2004, 2008, 
Hollenbach & Gorti 2009, hereafter GH04, GH08 and HG09 respectively) to compute gas emission 
line luminosities from the disk around TW Hya. In §2, we provide a short description of the available 
data and current understanding in each case. We briefly describe the modeling procedure and then 
discuss disk models that best match observational data in §3. We consider the implications of our 
results on the nature of the transition disk around TW Hya (§4), discuss the constraints set on disk 
mass in §5 and finally present our conclusions in §6. 



2. Gas Emission Lines: Observational data 

We summarize below the observed gas line emission at various wavelengths. For clarity of 
discussion, we do not present the data in chronological sequence of detection, but instead group 
observations into three wavelength regions, (i) far-infrared and sub-millimeter features, (ii) the mid- 
infrared (MIR) features, and, (iii) the optical and near-infared (NIR) features. Table[l]lists the line 
luminosities for the observed gas emission lines. We do not discuss the dust emission observations, 
but refer the reader to C02, Eisner et al. (2006), Ratzka et al. (2007) and references therein. 

Sub-millimeter and far infrared emission: Detection of CO 2-1 pure rotational emission from 
TW Hya was first reported in the seminal paper on gas disk lifetimes by Zuckerman et al. (1995). 
Several molecular species (CO, HCN, CN, HCO+, DCO+) have been detected in the TW Hya 
disk at millimeter and submillimeter wavelengths using single-dish telescopes (van Zadelhoff et al. 
2001; van Dishoeck et al. 2003; Thi et al. 2004). Wilner et al. (2003) presented the first spatially 
resolved observations of the disk in the HCO"*" J=l-0 line. This was followed by CO rotational line 
images in the 2-1, 3-2 and 6-5 transitions using the SMA by Qi et al. (2004, 2006). 

Thi et al. (2010), using the Herschel Space Observatory, recently reported the detection of the 
[OI]63^m line and set upper limits on the [OI]145/im and [CII]158/im lines. In conjunction with 
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^^CO measurements, they interpret the emission hne fluxes as indicative of a low gas/dust mass 
ratio in the disk (~ 2.6 — 26) and estimate the disk gas and dust masses as (0.5 — 5) x 1O~'^M0 and 
1.9 X 1O~'^M0 respectively. 

Mid-infrared emission: The [NeII]12.8/im forbidden line was first reported by Ratzka et al. 
(2007) from Spitzer IRS observations along with HI 6-5 and 7-6 emission lines. The [Nell] line was 
subsequently spectrally resolved using MICHELLE on Gemini N by Herczeg et al. (2007) and they 
found the line profile to be consistent with a single Gaussian centered at —2 it 3 km s~^ in the 
stellocentric frame with an intrinsic FWHM of ^ 21 km s~^. They concluded that the broad line 
width could result from turbulence in a warm disk atmosphere, Keplerian rotation of gas located 
~ 0.1 AU from the star, or a photoevaporative flow. These observations were followed by high- 
resolution spectra using VISIR on Melipal/VLT by Paccucci & Sterzik (2009) and more recently by 
Pascucci et al. (2011). They found that the [Nell] line proflle (intrinsic FWHM ~ 13km s~^) and 
the ~ 4 — 6 kms~^ blueshift with respect to the star are consistent with theoretical predictions of a 
photoevaporative flow driven by stellar extreme-ultraviolet (EUV) photons (Alexander 2006). The 
line profile may also be consistent with photoevaporation by X-rays and the production of [Nell] 
in predominantly neutral. X-ray heated layers (Ercolano &: Owen 2010). 

Najita et al. (2010) recently presented high quality Spitzer IRS spectra of TW Hya in the 
10 — 20/im region. The spectra show many emission features including H2 S(l) and S(2) pure 
rotational lines, many HI lines, [Nell], [Nelll], several transitions of OH and bands of HCO"*" and 
002- They find no strong lines from H2O, C2H2, or HON as seen in other classical T Tauri stars 
without dust holes (Carr &: Najita 2008, Salyk et al. 2008, Pascucci et al. 2009, Pontoppidan et 
al. 2010). They conclude that the [Nell] emission is consistent with origin in the X-ray heated disk 
surface. The detected OH emission arises from transitions originating in energy levels as high as 
23000 K above ground, and they interpret the emission as likely to arise from photodissociation of 
water in the disk since this process results in OH in very high rotational states (van Dishoeck & 
Dalgarno 1984). The relative strengths of the H I emission are found to be consistent with case B 
recombination in gas possibly associated with the accretion flow close to the star. 

Near-infrared, optical and UV emission: Emission from molecular hydrogen in excited vi- 
brational states has been detected from TW Hya. Weintraub et al. (2000) detected emission from 
H2 in the disk around TW Hya in the 1-0 S(l) line at 2.1218 fim through high-resolution spec- 
troscopy using CSHELL on NASA's Infrared Telescope Facility. Herczeg et al. (2002) detected H2 
UV emission (with FWHM~ 18km s~^) consistent with excitation caused by UV- induced fluores- 
cence from a warm surface layer ~ 2 — 3AU from the central star. Ingleby et al. (2009) identifled 
a feature in the STIS spectra at 1600 A which they attribute to electron impact H2 emission and 
derive a lower limit on the surface density of gas at ~ 1 AU to be ~ 5 x 10~^g cm~^. 

CO rovibrational emission from hot gas, presumably in the inner (r < 4AU) disk, has also been 
detected. Rettig et al. (2004) first detected CO from the disk, and concluded that the observed 
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emission was optically thin and consistent with a CO gas mass of ~ 6 x 10 g. The low emission 
temperature of ~ 430K was interpreted as due to a clearing of the inner disk out to a radial distance 
of ^ 0.5 AU. Salyk et al. (2007) expanded on these results with observations that cover a larger 
range of rotational energies. They estimate emission radii for the CO gas to be ~ 0.2 — 1 AU, and 
derive higher temperatures (~ 800K) based on their excitation diagram analysis. Pontopiddan et 
al. (2008) have recently presented velocity-resolved spectroastrometric imaging of the rovibrational 
lines (4.7/im) using the CRIRES high-resolution infrared spectrometer on the Very Large Telescope 
(VLT). Keplerian disk models are fit to the spectroastrometric position- velocity curves to derive 
geometrical parameters of the molecular gas and the CO gas emission from TW Hya is calculated 
to come from r ~ 0.1 — 1 AU. 

Optical lines from atomic oxygen have been detected from TW Hya. Emission from forbidden 
lines of [01] at 6300A and 5577A have been observed around many young stars with disks and 
outflows (e.g., Hartigan et al. 1995). Apart from a high velocity component (HVC, Av ~ 100 — 
200km s~^) that is attributed to shocks in the outflows, there is also a low velocity component 
(LVC, Av 10km s~^) that is believed to originate in the disk or from the base of a slow disk 
wind. The origin of this line has remained somewhat of a mystery, as physical conditions in disks 
are such that reproducing the observed high line luminosities is difficult (e.g.. Font et al. 2004, 
Hollenbach & Gorti 2009). The higher excitation [01] 5577 A hne often accompanies the LVC, but is 
never associated with the HVC. The optical forbidden lines often exhibit small blueshifts (typically 
~ —5 lb 2km s~^) (Hartigan et al. 1995). Ercolano &: Owen (2010) explain the LVC emission as 
arising from an X-ray heated, neutral photoevaporative flow that then gives rise to the observed 
blueshifts. However, the forbidden line emission in TW Hya is not blueshifted and does not appear 
to participate in a photoevaporating flow (Pascucci et al. 2011). TW Hya has emission lines (only 
in the low velocity component) at [OI]6300A and 5577A, with line luminosities estimated at 10~^ 
and 1.4 x 10~^ Lq respectively (Alencar & Batalha 2002, Pascucci et al. 2011). Pascucci et al. 
(2011) obtain an intrinsic FWHM of the [01] lines as ~ 10km s~^ and also report [SH] 4069A 
emission, estimated at 3.6 x 10~^Lq and non-detection of the [SH] 673lA line. As we show later, 
very high gas temperatures or densities are needed to excite the [01] 5577A line and reproduce 
the observed ratio of the 6300/5577 hues {Tgas > 30,000K for Ue < lO^cm'^, or, Tgas > 4000K 
for Ue ^ 10^ cm~^). Thus it is a challenge for the 5577A emission to arise thermally (however, 
also see Ercolano & Owen 2010). Here we account for the [01] emission and line ratio as a result 
of photodissociation of OH in the disk surface layers (e.g., Storzer &: Hollenbach 1998, Acke et al. 
2005). 



3. Gas Emission Lines: Modeling 

We use our recently developed theoretical disk models to model the gas emission from the disk 
around TW Hya. The stellar parameters are reasonably well determined and the dust distribution 
has been previously modeled, and the focus of this paper is therefore the gas disk modeling. We 
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vary disk parameters to find the model that best fits the observational data, with the principal 
variable being the surface density distribution. We find that two disk models with different inner 
disk masses (for r < 4AU), but with the same outer disk mass (r > 4AU), provide reasonable fits 
to the data. We describe these disk models in what follows, and discuss in detail the case with the 
lower inner disk mass as these results are a somewhat better match to the observed line emission. 
We call this our standard disk model (hereafter referred to as SML, 'L' for low inner disk mass). 
We also discuss the results from an alternate model with a higher mass in the inner disk and refer 
to that composite model as SMH. 



3.1. Model description and inputs 

Our theoretical disk models start from a prescribed disk surface density distribution and stellar 
spectrum and proceed to self-consistently calculate the disk structure from the gas temperature and 
density distribution. The gas disk models simultaneously solve for the thermal and chemical prop- 
erties of the gas and vertical hydrostatic equilibrium. Gas and dust temperatures are computed 
separately and can deviate significantly in the surface layers where most of the line emission origi- 
nates. The gas radiative transfer follows an escape probability formalism. Dust radiative transfer is 
done using a modified two-layer approach (e.g., Chiang &: Goldreich 1997, Dullemond, Dominik & 
Natta 2001, Rafikov & DeColle 2006). Our chemistry consists of ~ 85 species and ~ 600 reactions, 
includes photodissociation and photoionization by EUV, FUV and X-rays, and is focused towards 
accuracy in the disk surface layers. Our standard models do not include freeze-out of species on 
grain surfaces, but we test for the effects of freeze-out where relevant below. We also assume that 
H2 forms on grain surfaces in our disk models for all grain temperatures below ~ 400K (Cazaux & 
Tielens 2004). For further details, see GH04, GH08 and HG09. 

Stellar parameters and the incident radiation field are key inputs to our theoretical models, 
and are reasonably well determined for TW Hya. The source is at a distance of 51 pc (Mamajek 
2005) and the disk is almost face-on with a very low inclination angle (7° - outer disk, Qi et al. 
2004; 4° - inner disk, Pontoppidan et al. 2008). We adopt M = O.7M0, R=1R0 and T = 4200K for 
the stellar mass, radius, and effective temperature (Webb et al. 1999). We construct a composite 
spectrum for the star by adding observed X-ray and FUV spectra. We use the FUSE/IUE spectrum 
from the NASA HEASARC archive (Valenti et al. 2003, Herczeg et al. 2002). The FUV region is 
split into several energy bins for accurately computing gas opacity as described in GH04. Herczeg 
et al. (2002) derive a reconstructed Lyman a flux ~ 30% higher than the observed flux and we 
use this value for the Lyman a band of our FUV spectrum. The total FUV (6 — 13.6eV) flux is 
approximately 3 x 10^^ erg s~^. The X-ray spectrum is from the XMM-Newton archive (Robrade 
& Schmitt 2006) and extends from 0.3 — 5keV, and with a total X-ray luminosity of 1.3 x lO^'' erg 
s~^. We extrapolate the observed spectrum to lower and higher energies (0.1 — lOkeV) and this 
results in a total X-ray luminosity of 1.6 x 10^° erg s~^. Although the star might also have an EUV 
component (13.6 — lOOeV), perhaps indicated by the observed photoevaporative flow (Alexander 
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2008, Pascucci & Sterzik 2009), its strength can only be crudely estimated by indirect arguments 
(Alexander et al. 2005) and we do not include this uncertain contribution to the high energy flux. 
We discuss below, where relevant, the effect on line emission fluxes if an EUV flux were to be added. 
Figure [J shows the composite spectrum used as model input and the vertical dashed lines in the 
figure indicate the neglect of an EUV radiation field. 

Two different dust disk models have been proposed for TW Hya, by Calvet et al. (2002) 
and Ratzka et al. (2007). C02 developed detailed models of the TW Hya dust spectral energy 
distribution (SED) to constrain many aspects of the disk structure, and these models have been 
validated by subsequent observational studies (e.g., Qi et al. 2004, Eisner et al. 2006, Hughes et 
al. 2007). The flux deficit at < 20;um and a sharp increase or excess at ~ 20 — 60/xm led them to 
propose the presence of an inner hole in the disk at r < 4 AU. The deficit is ascribed to a lack of 
dust in the hole, and the excess arises from direct illumination of the rim of the optically thick dust 
disk beyond 4 AU. There is also a small mass of dust, 0.5 lunar masses, in the inner cavity which 
gives rise to the observed spectral feature at 10/im. Ratzka et al. (2007) situate the inner edge of 
the outer optically thick dust closer to the star at ~ 0.5 — 0.8 AU. They model both the SED and 
the 10/im visibilities in detail by including sedimentation in the disk, where the dust is reduced in 
the upper atmosphere due to settling processes. Our dust models at present do not include settling. 
We defer gas modeling of the disk using the Ratzka et al. (2007) dust distribution to future work, 
and keep our dust model relatively uncomplicated at present. 

We avail ourselves of the dust models of C02. For simplicity, we adopt many features of 
this dust model, specifically the dust size and spatial distribution in our standard dust model. 
We consider their favored model with a maximum grain size amaa:=lcm and disk mass O.OGMq, 
which corresponds to a total dust mass of 7.8 x 10~^Mq for their assumed gas/dust ratio. As we 
ignore dust particles with sizes larger than 1mm which are not very relevant for gas physics, our 
model disk "dust" {amax = 1mm) mass is reduced accordingly to 2.4 x lO~^M0(for a grain size 
distribution n(a) oc a~^'^, where a is the grain size). Our chemical composition of dust for the 
outer disk is simpler than the detailed models of C02. We use a mixture of silicates and graphite 
for the optically thick disk, and verify the assumption by the SED fit from the resulting model dust 
continuum emission. We note the disk structure (especially in the regions where gas emission lines 
originate) is eventually determined by the gas temperature structure, although the dust distribution 
is sometimes important for heating and cooling of gas. We use pyroxene grains for the inner dust 
hole as suggested by C02 to better fit the continuum SED and the lO/xm feature. As we show later, 
the amount of dust present in this inner region is in any case too low to significantly affect the 
gas thermal balance or chemical calculations. The details of the adopted dust model are listed in 
Table M 

Polycyclic aromatic hydrocarbons (PAHs) can efficiently heat gas via the grain photoelectric 
heating mechanism (Bakes &: Tielens 1994) and gas heating depends on the assumed PAH abun- 
dance. Disks around low mass stars are believed to be deficient in PAHs by a factor of ~ 10 — 100 
(Geers et al. 2006). We typically scale the PAH abundance with the decrease in cross-sectional area 
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of the dust per H atom compared to the interstehar value. We adopt a PAH abundance of ~ 10~ 
per H in the ISM (e.g.,Habart et al. 2003). Therefore, the abundance of PAHs in the optically thick 
disk (r > 4AU) around TW Hya is ~ 10~^. Such a depleted abundance means that X-ray heating 
usually dominates FUV-induced PAH heating of gas in the outer disk. The abundance of small 
grains (a < 0.1/um) in the inner disk (r < 4AU) is either negligible (C02) or very small (Eisner et 
al. 2006). It is therefore likely that the abundance of very small dust like the PAHs is low in the 
inner disk. Such low PAH abundances do not influence the heating in the inner disk and the PAH 
abundance in the inner disk is therefore set to zero. 

We do not consider turbulence in the disk beyond 4 AU and all cooling line widths are thermal 
in the optically thick outer disk. This assumption is supported at large radii by sub-millimeter 
observations and modeling of molecular emission from TW Hya by Qi et al. (2006) and Hughes 
et al. (2010). Hughes et al. (2010) infer a small turbulent linewidth that is approximately a 
tenth of the thermal linewidth in the outer disk. The situation may be different in the inner disk, 
r < 1 — 2AU. Salyk et al. (2007) assume higher turbulent linewidths {vturb ~ 0.05vk ~ 4km s~^, 
where the Keplerian velocity vk ~ 80km at 0.1 AU) in their analysis of the rovibrational CO 
lines from the disk. COrovib emission arises from hot gas, and the thermal speed of CO for ~ 2000K 
gas is ~ lkms~^. The expected thermal linewidth is therefore ~ 2km s~^. Observed line widths 
are typically ~ 15kms~^ (Salyk et al. 2007), with presumably a large Keplerian contribution for 
gas located at small radii (~ O.lAU) in a face-on disk with a small inclination (4°, Pontopiddan et 
al. 2008). We follow Salyk et al. (2007) in assuming Vturb = O-O^vk in the inner disk (r < 4AU) 
regions. 

We also add viscous heating as a mechanism to heat gas in the inner disk. The dust-depleted 
inner cavity might still contain dense gas and viscous heating may be an important contributor 
to heating the dust-poor gas, especially in the shielded interior near the midplane. We follow the 
prescription given in Glassgold et al. (2004) (Equation 12 of that paper) for the viscous heating 
term. The viscosity parameter a is unknown, but is set in our models by assuming steady state 
accretion onto the star at the observed rate of ~ 10~^ M0yr~^ and from the solutions obtained for 
the density and temperature of the gas. 

The surface density distribution of gas in the disk is the primary input to our disk models. 
We begin with the dust distribution of C02 which has a rapid increase in surface density over 
~ 3.5 — 4 AU. We then vary the gas surface density distribution, S(r), so that the resulting gas 
emission lines provide a good match to observational data. The resulting profile for the SML disk 
is shown in Figure [2j The SML disk has a gas mass ^ 10~^ Mq within 3.5AU, and O.O6M0 for 
3.5AU < r < 200AU. Although the dust surface density drops by ~ 10'^ inward of r 4AU, the gas 
surface density drops only by ~ 100 for the SML disk. The alternate model SMH has a gas mass 
~ 10"^ M0 within 3.5AU, and in this case the gas surface density only drops by ~ 10. The outer 
disk surface density distribution is the same for both SML and SMH cases. We discuss later in §3.3 
how deviations from the models SML and SMH worsen the agreement between model results and 
observational data. 
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For descriptive purposes, we divide the disk into three distinct regions, (i) the inner disk 
(r < 3.5AU), (ii) the mid-disk region (3.5AU< r < 20AU) and the (iii) outer disk (r > 20AU). 
Our modehng procedure first constrains the gas in the inner disk region, as this gas provides 
opacity to stehar radiation, shields gas in the partiaUy exposed rim region at 3.5 — 4 AU, and 
affects the structure of the outer disk. The inner region of the SML disk is not only depleted in 
dust but also in gas, by > 2 orders of magnitude compared to that expected if the outer optically 
thick disk was smoothly extrapolated radially inwards. A power law is assumed for the inner disk 
surface density distribution, S(r) oc r~^; 0.06AU< r < 3.5AU. We place the inner disk edge at 
0.06AU, at the location of the magnetospheric truncation radius (Eisner et al. 2006). As there 
is continued accretion onto the star, presumably by funnels along magnetic field lines, this is a 
logical choice for this parameter. Between 3.5 — 4 AU, S(r) rises sharply and these radii constitute 
the "rim" region. The functional form for S(r) is modeled to fit the appearance of the rim of a 
photoevaporating disk (Gorti et al. 2009) and is given by logi;(r) = 18 — exp(ll/r^[/) g cm^^; 
3.5 < rAu < 4. Here, rAu = t/1AU. The dust disk turns optically thick over this region. Beyond 
4AU, the surface density is assumed to follow a self-similar density profile, as expected from a 
photoevaporating and viscously evolving disk (e.,g., Hughes et al. 2007, Gorti et al. 2009). For 
r > 4AU, S(r) = 500r;^[^^ exp(-ri^/100) g cm"^. 

3.2. Model Results: Disk Structure 

We briefly discuss the disk density and temperature structure, the dominant heating and 
cooling mechanisms at three representative disk radial positions, and then describe the line emission 
calculated from the disk models. 

The density and temperature distribution as a function of spatial location in the SML disk 
is shown in Figure [3l Gas temperatures in the surface regions where emission lines originate are 
typically ^ 2000 — 200K in the inner and mid-disk regions, r < 30 AU. The near absence of dust 
in the inner disk and the sudden rise in dust surface density at the rim causes an increase in gas 
heating and results in a small vertically extended region at this radius ~ 4AU. However, gas in 
the inner disk provides opacity that shields the midplane regions at ~ 4 AU, and there is only a 
slight "puffing-up" of the inner rim. (We find that for the upper layers of the outer disk where the 
gas temperature decouples from the dust temperature, the column density of gas to the star at a 
given angle (height of the rim) changes by less than 10% as the mass of the inner disk changes from 
1O-^M0 to IO-'^Mq.) As seen in Fig. [3] (see Ay = 1 surface), the disk structure is overall quite 
flat, with only modest flaring of the outer disk. Gas at depths where Ay ~ 0.1 — 1 to the star, and 
at radii r > 50 AU, only attains temperatures of ~ 50 — 200K and the density in these warm layers 
is low, n ^ 10^"'^ cm~^. 

Disk gas temperature and density is shown as a function of height at typical positions in the 
three different disk regions, inner disk at r = 0.1 AU, mid-disk at r = 6AU, and the outer disk at 
r = 30AU in Figures HIS] and El 
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In the inner disk (Fig. heating is mainly by X-rays all through the disk, and in the midplane 
FUV is significant only via chemical heating (primarily due to photodissociation of water and 
ionization of carbon). However, to attain the high surface temperatures achieved in this region, 
FUV photons are essential to photodissociate molecules capable of efficiently cooling the surface gas. 
The gas temperature increases when H2 photodissociates and molecular (H2O and OH, formed via 
H2) cooling decreases, and also when CO photodissociates and can no longer cool the gas efficiently. 
In the inner disk, the transition from C to CO occurs higher up than the transition from H to H2, 
as is typical in low dust regions (e.g., Glassgold et al. 1997). The inversion in height can be 
seen in Fig. H] where the H2 and CO photodissociation fronts are marked. Rovibrational emission 
of CO {COrovib) dominates the cooling with some additional cooling by H2O and OH, until CO 
photodissociates. At higher z, OI6300A and Lyman a are the main coolants. 

Fig. [5] shows the disk vertical structure in the mid-disk at r = 6AU. In these regions where 
dust is abundant, dust collisions keep gas at the ambient dust temperature in the optically thick 
midplane region. Near the optically thin surface layer and near the H2 photodissociation front, 
heat released from the formation of H2 is significant. Chemical heating due to photoreactions 
(FUV photodissociation of water) and FUV-induced grain photoelectric heating by PAHs is of 
significance over a limited column (Ay ~ 0.1 — 1). At the surface. X-rays dominate the heating. 
Cooling is mainly due to CO and H2 in the molecular layer, and OI6300A, [OI]63/im, [Nell] and 
[Aril] fine structure lines cool the surface gas. 

The outer disk vertical structure at r = 30AU is shown in Fig. O Dust is thermally coupled 
to the gas by collisions in the midplane where gas and dust temperatures are equal. Near the 
Ay = 1 layer, heating by PAHs, photodissociation reactions and formation and dissociation of H2 
are important. X-rays dominate the heating at vertical column densities of < lO^^cm"^. FUV 
photoionization of carbon contributes to the heating (~ 15%) in the region above the CO photodis- 
sociation layer. Main coolants in the outer disk are CO in the molecular region and [OI]63;Lim fine 
structure emission at the surface. 

3.3. Model Results: Line Emission 

We adopt the following procedure to determine the surface density distribution that is described 
in §3.1. In order to obtain a gas disk surface density model that best matches observed gas line 
emission, we first constrain the surface density distribution in the inner disk, i.e. within 4 AU. 
This is because gas in the inner disk can shield the outer disk regions, affecting the density and 
temperature structure of the disk beyond 4 AU, where the dust disk is optically thick. There 
clearly is some gas in the inner hole, based on the observed accretion onto the central star, and 
on observations of CO rovibrational lines and H2 ultraviolet lines. We discuss these emission line 
constraints and the undetected water lines below, and their implications for the gas mass in the 
inner hole. 
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Gas Opacity Although the dust is opticahy thin to stellar radiation, the column density of gas 
may be sufficient to provide significant optical depth in the inner disk. In fact, the gas may be 
optically thick to its own infrared emission. We use our disk chemical models to obtain chemical 
abundances for gas in the inner disk and follow Tsuji(1966) in calculating the monochromatic 
opacity due to each gas species. At the likely pressures and temperatures of our inner disks, H~ 
and H2 bound-free and free-free, CO and H2O all can provide significant opacity. H~ and 
bound-free and free-free continuum opacity can be significant at wavelengths shorter than about 4 
Hm, and produce near IR radiation from the inner edge of the disk at ~ 0.06 AU. The ground state 
CO vibrational band occurs near 4.7 fiui and the blend of CO rovibrational lines provide opacity 
near 5 //m. Beyond about 6 /xm, the blend of pure rotational H2O lines dominate the opacity. 

Later on in this subsection, we discuss CO and H2O lines in some detail. Here, we make more 
general arguments based on the expected continuum emission from gas in the disk. We use the 
opacities from Tsuji (1966) to show that if the surface density of the outer disk is extrapolated 
inwards to produce a ~ 10-3 Mq inner disk (we call this the "full disk" or the "undepleted disk"), 
observable emission would result in the wavelength region 5 — 15 /im due to the CO and H2O 
emission. We find, by a similar analysis, that the H" and bound-free and free-free produce 
excess emission at shorter wavelengths (~ 2 — 3/im) where the stellar emission is significant. We 
use the specific opacity (per molecule) for each opacity contributor from Tsuji(1966) and calculate 
the opacity at each radial annulus by multiplying this with the vertical column density as obtained 
from our model disk chemical calculations. For illustration, we calculate the flux at one wavelength 
and compare that with the observed emission. We choose a wavelength of 8/xm where the stellar 
contribution is negligible, there is no contamination from the silicate feature and where there is a 
significant dip in the continuum emission (see C02). At each radial annulus, the photosphere (at 
8//m) is defined to lie at the height (z) above the midplane where the total optical depth at 8/im 
becomes of order unity. Near the inner edge of the disk for 0.06AU < t < 0.08AU, H" and 
free- free are the main opacity contributors. This opacity drops rapidly with decreasing electron 
density further out in the disk. At 0.08AU ^ r < 0.2AU, H2O lines dominate the 8^m opacity. 
Having defined the photosphere of the disk from the gas opacities, we estimate the 8fj,m fiux from 
this surface, where the gas temperature ranges from ~ 300 — 2000K. The disk becomes optically 
thin to the midplane (A ~ 8/zm) at ~ 0.2 — 0.3AU. We obtain an 8//m fiux of ~ IJy, a factor 
of ~ 4 higher than observed. A full gas disk is therefore incompatible with the spectral energy 
distribution. The H2O lines at 8/im therefore yield an estimate of the upper limit on gas mass as 
2.5 xlO"^ Mq. We discuss later our model results for H2O line emission and compare them to 
Spitzer IRS upper limits to the line fluxes to obtain a better limit on the inner disk mass. 

A similar opacity analysis at a wavelength of 5^m suggests that CO rovibrational emission in 
excess of what is observed will be produced by a full gas disk. We use our model calculations and 
compare them with observed rovibrational CO line strengths near 5/im to estimate the actual mass 
of the inner disk to be ~ 10"'' - 10"^ Mq. 

In what follows, rather than using an opacity argument, we use our disk model results on 
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emergent line fluxes from the disk surface to infer the disk mass. We describe the emitting regions 
for each gas emission line diagnostic and also discuss where applicable the constraints set in each 
case. The outer disk gas is attenuated by the gas in the inner disk as described earlier. We calculate 
the outer disk structure using the surface density distribution of the inner disk in the SML model 
and use the same results for the outer disk component of the SMH model as well, as the attenuation 
columns of the gas provided at 4 AU are similar in both these cases (within '^10%). As our disk 
models are computationally intensive, we do not calculate disk models for the outer disk in the 
SMH case separately. In the rest of the paper, the emission in the mid and outer disk regions refer 
to that calculated using gas optical depths and absorption columns for the 10~^Mq inner disk mass 
case (SML). Figure [7] shows the radial extent of each of the gas emission lines discussed below. The 
cumulative line luminosity, L{r), integrated up to the radius r is shown for each emission line and 
is normalized to the total calculated value (from inner, mid and outer disk regions). Also marked 
are the 10% and 90% luminosity levels. We quote the corresponding radii below for each emission 
line and this defines the radial extent of the emission. 

CO rovibrational emission COrovib emission typically arises from a height in the inner disk 
where the CO is just beginning to photodissociate. Gas temperatures in the disk generally increase 
with height and gas temperatures in the emitting layers are high {Tgas ^ 500K). In dust-poor gas, 
the CO and II2 photodissociation layers are inverted, as first reported by Glassgold et al. (2004), 
that is, the carbon becomes CO at larger z than the H becomes II2. We similarly find that the 
II2 is photodissociated to deeper layers and that CO is co-existent with large columns of atomic 
hydrogen, creating conditions favorable for exciting COrovib emission (Glassgold et al. 2004). The 
primary formation route to CO is the oxygenation of atomic carbon (C-I-OH). In the layers shielded 
from FUV photodissociation CO is destroyed by collisions with helium ions. 

CO vibrational lines are an important constraint for setting the gas mass in the inner disk. 
A high disk mass results in high total CO vibrational emission if the lines are optically thick, 
as the luminosity in these lines increases with temperature which increases (albeit weakly) with 
increasing gas mass. The total luminosity of the observed lines (select P and R branch transitions) 
is constrained by observations to be < IO^^Lq (Salyk et al. 2007). We find that inner disk masses 
of 10"'*, 10^^ and 10"^ Mq all give total CO luminosities within a factor of ~ 2 of the observations 
and are all hence viable choices. The CO emission is found to be mainly a surface effect dominated 
by hot gas at the surface of the disk, resulting in approximately the same total CO luminosity 
for all three disk masses. To distinguish between these cases, the line luminosity of each of the 
CO rovibrational lines needs to be compared with data. We estimate the CO rovibrational line 
luminosities of the v = 1 — P(J) series as follows. We adopt the analysis of Scoville et al. (1980) 
and perform an approximate non-LTE calculation for the populations of the CO rovibrational 
levels. We only consider the v = and v = 1 states as higher vibrational states are unlikely to be 
significantly occupied at typical disk temperatures. We assume that collisions result in a change 
in the vibrational state, and vibrational level populations are determined from detailed balance 
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assuming an "effective" Einstein A- value from u = — 1 as described in Najita et al. (1996). We 
then assume that the rotational levels within each v state are thermally populated (Scoville et al. 
1980). Collisional rates with e~, H and H2 are taken from Najita et al. (1996). Oscillator strengths 
and transition frequencies are from the compilations of Kirby-Docken & Liu (1978). The optical 
depth for each transition is calculated as detailed in HoUenbach &l McKee (1979) from the column 
density of CO molecules in a given rovibrational state, and an escape probability formalism is used 
for line transfer. 

Salyk et al. (2007) detect 15 transitions covering rovibrational energies E^/k from 3000K- 
6000K, where is the energy of the excited upper level. We also compare our results with earlier 
CO rovibrational detections by Rettig et al. (2004) in complementary transitions. There was a 
typographical error in the Rettig et al. (2004) table of line fluxes, and the actual fluxes are a factor 
of 10 below what has been reported. We also do not include the P (2 1-25) transitions of that paper, 
as these data are prone to errors due to photospheric and telluric corrections (Sean Brittain, private 
communication) . We compare our model results with the data in Figure El We find that a gas 
depletion from the full disk of 0.1 — 0.01 corresponding to a gas mass of 10"^ - lO-^M© matches the 
observed line emission fairly well, with the lower mass inner disk model (SML) being a better fit to 
the CO data. A factor of 10 lower in mass results in CO emission from warmer gas that produces 
excess emission at high J levels. All of the CO gas emission comes from r < 0.5AU (due to rapidly 
decreasing gas temperatures beyond), in agreement with both the analysis of Salyk et al. (2007) 
and the spectroastrometric imaging study by Pontopiddan et al. (2008). Our gas mass estimate 
for the SML disk is consistent with the surface density estimate of Salyk et al.(2007) of 1 g cm~^ 
at r ~ 1 AU where they used a slab model, although in our models the emission is mainly from hot 
gas in the disk surface layers. Salyk et al. (2007) also assume a H2/CO conversion ratio of 5000 
as observed towards a massive protostar, and caution that their derived mass estimate is likely 
be a lower limit. In our disk models, the hydrogen is predominantly atomic in the regions where 
COrovib emission originates (also see Glassgold et al. 2004). We explicitly calculate disk chemical 
abundances and find that the H/CO ratio in the hot inner disk regions where COmb emission arises 
is ~7000. The gas/dust mass ratio in the inner disk is ~ 500 — 5000 for the dust mass of the Calvet 
et al. (2002) disk, indicating a relative depletion of the dust with respect to the gas at r < 4AU. 

H2O rotational emission Our gas opacity analysis above suggests that a full undepleted disk 
would produce excess H2O rotational emission in lines whose wavelengths range from 6 to 15 //m 
compared with observations. As we specifically calculate water abundances and the disk density 
and temperature structure, we can compare the model disk H2O rotational line fluxes to the 
corresponding upper limits from Spitzer IRS in the 10 — 19/im wavelength region (Najita et al. 
2010). No water emission has been detected and line luminosities are therefore expected to be 
lower than ~ 10~^Lq. We find that for the full undepleted inner disk with a mass of 1O~^M0, 
many rotational lines in the 10 — 19/im region are above the detection limit. Some of the strongest 
transitions are 7qi — 6i6,762 — 615,854 — 707,853 — 726 and 844 — 7i7, all with luminosities in the 
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range of 1 — 3 x 10~^Lq, compared with the upper hmits set by Spitzer observations of ~ 1O~^L0. 
From the water rotational hne emission constraint alone, we can set a limit on the inner disk mass 
as ~ 1O~^M0, at which the water emission is barely detectable. This limit is consistent with the 
upper limit set by CO emission in the preceding discussion, where we fix the inner disk mass to be 
in the range 10~^ — 10~^Mq based on CO rovibrational emission. 

H2 fluorescent emission We do not model H2 UV fluorescence but draw on the results of 
Herczeg et al. (2004) who conduct a very detailed analysis and model the observed emission. Our 
SML disk corresponds to their thick disk model, and most of the emission in this case arises from the 
illuminated edge of the disk at r = 0.06AU. Herczeg et al. (2004) estimate a hot {Tgas ~ 2500K), 
H2 mass of ~ lO^^g from their analysis. We calculate the mass of H2 gas at these temperatures in 
our inner disk to be ~ 3 x lO^^g, in reasonable agreement with the Herczeg et al. (2004) results. 
All of this gas lies very close to the inner edge at 0.06AU. At this radius the Keplerian speed 
of the gas from the almost face-on disk {i = 4°) is about 7km s~^, the thermal speed is ~ 4km 
s~^, bringing the expected FWHM to ~ 16km s~^ very close to the observed FWHM of 18km s~^ 
(Herczeg et al. 2002). We note that the inner disk structure near the magnetospheric truncation 
radius is likely to be quite complicated compared to our simple power law assumption of the surface 
density, and our comparison with the Herczeg et al. model is only a consistency check. 

OH lines 16 rotational transitions of OH, mainly in the ground vibrational state, were reported 
by Najita et al. (2010) from their Spitzer observations. The data indicate high excitation states 
{3upper ^ 15-30), with equivalent energies E^/k up to 23,000K above ground. OH mid-infrared 
thermal emission from our model is mainly from the hot, dense (n ~ lO'^^ — lO'^^cm"^, T^^s ~ 
700 — lOOOK) gas very close to the central star at r < 0.1 AU. These lines are, however, weak in the 
10 — 19/im region and below the sensitivity of the Spitzer IRS instrument. OH thermal emission is 
also produced from gas in the mid-disk and outer disk regions, but this gas is typically at ~ 300K 
and incapable of reproducing the observed high excitation lines. As suggested by Najita et al. 
(2010) for TW Hya, UV photodissociation of water is a possible source of the OH emission (see 
also Tappe et al. 2008 who suggest this mechanism for HH214). 

From our disk chemical network solution, we can estimate the total luminosity in the OH lines 
produced by the photodissociation of water in the disk. We calculate the total OH line luminosity 
as follows. For simplicity, in this discussion we assume that Lyman a photons {E = 10.2eV) are 
mainly responsible for the photodissociation of water molecules - a reasonable assumption because 
Lyman a contributes to a significant fraction of the total FUV luminosity (e.g., Herczeg et al. 
2002, Bergin et al. 2003). Harich et al. (2000) studied the photodissociation of water by Lyman 
a photons and the photodissociation energy is determined to be 5.1eV. The excess energy for a 
Lyman a photodissociation is therefore 5.1eV. About 66% of the photodissociations of water result 
in OH and H, while the remaining 34% result in O and H atoms. Harich et al (2000) further find 
that most of the OH products are extremely rotationally excited with a peak at a rotational level 
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J = 45 {Eu/k ~ 45,000K) and that almost 75% of the available energy goes into pure rotational 
excitation. This translates to an energy of ~ 4eV per photodissociation (of H2O to OH) that is 
then eventually radiated away in a rotational cascade from the high J states (also see Mordaunt et 
al. 1984, van Harrevelt & van Hemert 2000). We calculate the energy in the rotational cascade in 
our model disk by counting the total number of water photodissociations by photons with energy 
equal to or greater than Lyman a that lead to OH and multiplying this by the energy available 
for excitation (4eV). If Np£,{H20) is the rate of the total number of photodissociations of water in 
the disk by Lyman a (and more energetic FUV) photons, then the luminosity in the OH cascade 
LoH,PD is the product 0.66 x Np£){H20) x 4eV. Equating this to the observed OH luminosity 
then requires that there are ~ lO'*'^ photodissociations per second, or that only ~ 1% of the stellar 
Lyman a luminosity of 10^^ photons s~^ be intercepted by water in the disk (also see Bethell & 
Bergin 2009). If the solid angle subtended by the disk is ~ 0.1, the desired Lyman a flux absorption 
by water is feasible if water abundances are high enough (at the 10% level) to compete with the 
other dominant sources of opacity such as dust, photoionization ionization of Mg, Fe, and Si, and 
photodissociation of O2 and OH. Therefore, photodissociation of water can plausibly account for 
the observed line luminosities. 

In the above approximate analytic calculation, we simply used the Lyman alpha flux. However, 
in the model we calculate the OH luminosity due to water photodissociation using photodissociation 
rates obtained at each spatial grid cell. We use the full FUV spectrum, including Lyman alpha, 
but only count those photodissociations where the FUV photons are energetic enough to lead to 
OH in highly excited rotational states as described above. The OH luminosity in the SML disk is 
4 X 10~^Lq, lower than the total luminosity in the observed lines (upper J levels ranging from ~ 30 
to ~ 15), 6 X 10~^Lq. We note that the cascade must also produce lines that are not observed 
in the Spitzer IRS band. The energy of a rotational state oc J( J + 1) and the cascade originates 
at J ~ 45. The unobserved transitions from J ~ 45 to J ~ 30 and J~15toJ = therefore 
have a total energy approximately 1.5 times that in the transitions that are observed. The unseen 
luminosity must account for another ~ 9 x W~^Lq, bringing the "observed" value to 1.5 x IO^^Lq. 
Our SML disk results for OH produced by H2O photodissociation are discrepant by a factor of ~ 4. 
Of this emission, 30% arises from photodissociation of water in the inner disk and the rest comes 
from the extended middle and outer disk regions to radii r < 40AU (Fig. [7]). In the inner disk, the 
OH is from the edge of the disk near the star where the number densities are ~ cm 
and gas temperatures are ~ 500 — lOOOK. The main formation route to water is H2 + OH and it 
is destroyed mainly by FUV photodissociation. The water abundance in these regions is however 
low, with typical vertical column densities < 10^^ — 10^^ cm~^ in the photodissociation layer. In 
the mid-disk and outer disk regions, OH emission is dominated by the rim at 4 AU and beyond 
the rim is primarily produced from regions where the column density to the star ~ 1 — 5 x 10^^ 
cm~^, densities are 10^ — 10^*^ cm^'^ and Tgas ~ 200 — 600K. In the mid-disk regions OH-I-H2 
is the dominant formation route to water, whereas in the cooler outer disk regions (r > 15 AU), 
HsO"*" recombination leads to the formation of water molecules. Water abundances range from 
X{H20) ^ 10"^ to 10~^ in these photodissociation layers, and typical water columns to the star 
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are ~ 10^^ cm~^. The latter is consistent with our earlier analysis, which suggested that if the water 
opacity were ~ 10% of the total opacity at FUV wavelengths, we would obtain OH luminosities 
close to that observed. A column density of ~ 10^^ cm~^ of water provides optical depth of unity 
at the FUV wavelengths needed to photodissociate water. 

If the OH is mainly prompt emission following photodissociation of water, the SML disk 
appears to produce less water by a factor of ~ 4 than required by the observations. We consider 
a few mechanisms by which the water production rate may be enhanced to match observations. 
Raising the gas temperature and hence increasing the water formation rate by the endothermic 
H2+OH reaction is one possibility and this could be achieved by increasing the unknown PAH 
abundance and FUV grain photoelectric heating in the disk. However, as H2 pure rotational 
emission originates at the same spatial location where there is warm molecular gas, any increase in 
temperature overproduces the H2 S(l) line flux. We discuss the H2 S(l) and S(2) line fluxes below. 

Another possibility is that the SMH disk is the better solution for the inner disk mass. This disk 
is a reasonable match to the CO lines (Fig. [8|) and does not violate other observational constraints 
such as the lack of mid-IR water emission. Higher densities in the inner disk result in higher water 
abundances and the prompt OH emission increases by a factor of '--^ 4 in the inner disk, bringing 
the total OH luminosity to ~ 7.8 x 10^^ Lq, and only a factor of two below the observed value (see 
Tabled]). This model slightly overestimates the OI6300A and H2 S(2) line luminosities, but may 
nevertheless be a feasible model for the TWHya disk. The SMH disk also has the advantage that it 
better fits the [OI]6300A linewidth and provides a more typical value for the turbulence parameter 
a in the inner disk (see Table [6]) . 

Water formation on dust grains is another mechanism that could increase water in the disk. 
Our disk models do not include freeze-out and desorption processes to properly consider grain- 
surface reactions. We estimate an upper limit to the effect of water formation on dust in a simple 
manner as follows. We assume that oxygen atoms stick on grains, are instantly hydrogenated and 
released into the gas phase as water, at a rate equal to the collision rate of O with grains. We 
applied dust formation of gas phase water only in regions where the dust is cold enough to prevent 
thermal evaporation of O or OH before H2O is formed on the grain surface, Tg^ > lOOK. We use 
this relatively high grain temperature for thermal desorption so that we obtain an upper limit on 
the production of water on grains. We also require that Ay to the star is less than ~ 3 so that 
photodesorption can clear the water ice from the grain surface (Hollenbach et al 2009). Glassgold 
et al. (2009) considered a similar mechanism to increase water production in disks but found that 
it was not significant in increasing column densities of warm, observable water. However, the gas 
temperature is not relevant for our calculation as the OH lines are a result of photodissociation 
and not thermal emission. Water formation of grains in the SML disk (with IO'^Mq within 4AU) 
leads to a higher OH luminosity, 7 x 10~^Lq, but still a factor of ~ 2 below what is observed. 
We emphasize that this 7 x 10"^ Lq is an upper limit, given the optimistic assumptions on water 
formation on grains. 
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Radial transport of solids or water-bearing ices in the disk, a process not considered by us, has 
been suggested as a possible agent causing local enhancements and depletion of water abundances 
in the disk (Ciesla & Cuzzi 2006). Water ice from the midplane may also be transported to 
greater heights where it is subject to thermal desorption or photodesorption, leading to higher 
water abundances in the disk surface where it can be photodissociated. These effects, while likely 
to operate in disks, require time-dependent models that consider radial and vertical transport and 
are beyond the scope of this paper. However, one simple calculation is illuminating. If water is 
transported as water ice radially inwards to inner regions where it is photodissociated once, the 
mass influx needed is lO^'' H2O molecules per second. Assuming a H/H2O ratio of 2000 (if every 
O atom is in H2O ice) this translates to 2 x 10^^ H atoms per second or 5 X 1O-^M0 yr ^ . Since 
TW Hya has an age of lOMyr, and with a constant influx, this requires an improbable initial mass 
reservoir of 5Mq in the outer disk and also implies the current disk will only last 10^ years. This is 
unlikely. On the other hand, if turbulent mixing were to act vertically, then this has the advantage 
that the water can be reformed at the midplane and brought to the surface to be photodissociated 
repeatedly. For our adopted gas mass of O.OGMq, we therefore require that this mass be brought to 
the surface of the disk where it is photodissociated every 10^ years. At a typical radius, r = lOAU, 
the distance from midplane to the Ay = 1 surface layer where water is photodissociated is ~ lAU. 
Thus, the average vertical transport speed is only 1AU/I0^yrs= 5 cm s~^. This suggests that 
vertical mixing may supply the needed water. 

We have discovered a new mechanism for the rotational excitation of OH that may be the 
most promising, but difficult to quantify precisely at this time. Here, we do not need a greater 
production of gas phase water. The process is initiated by the photodissociation of OH, which leads 
to atomic oxygen in an excited electronic state (O^D) approximately 50% of the time (van Dishoeck 
&: Dalgarno 1984). The state of oxygen, due to its empty 2p orbital, is more electrophilic than the 
triplet ground state and is highly reactive to readily undergo bond-forming addition reactions. 
0(^D) + H2 re-forms OH very efficiently (with a rate coefficient 70 = 3.0 x 10^^*^ cm'^ s~^; NIST 
Chemical Kinetics database) leading to an effective reduction in the overall destruction of OH. 
This route has no thermal barrier as in the reaction of the ground state 0(^P) atoms with H2 
and is especially important at high gas densities when the rate of formation of OH by this route 
is faster than the radiative decay of the O^D atom to the ground state (with A ~ 8 x 10~'^s~^). 
The fraction fp of O^D that reacts with H2 before radiatively decaying to the ground state is 
given by fo = 1/(1 -|- {ricrit / ^{112))) , where the critical density ricrit ~ 7o/^ = 3 x 10''cm~'^. An 
important consequence of the 0(^D) -|- H2 reaction may be that this reaction produces OH in a 
rotationally excited state (with an efficiency eoH ~ 0.2 at J~20-30 states, Lin &; Guo 2008) which 
could then cascade to lower J to produce the observed MIR emission by Spitzer. In regions with 
n[H2) > ricrit, photodissociation of OH itself (rather than H2O as considered earlier) may result 
in the re-formation of rotationally excited OH. We estimate the contribution of this mechanism to 
the OH line emission as follows. 

Nqh photodissociations per second of OH will produce ~ O-beonfo^OH OH molecules in the 
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Jr^ 20 — 30 state per second. For our SML disk, we estimate that this mechanism could result in 
~ 5 X 10~^Lq (nearly all from the inner r < lAU region), which recovers most of the observed 
OH emission which the model was deficient in when we considered the photodissociation of water 
alone. We will explore this mechanism for producing OH emission in greater detail in a future 
study where we win treat O^D as a separate species for better accuracy, and treat the efficiencies 
of ^oh{v-, J) individually. 

We conclude with a footnote to the above discussion. Mandell et al. (2008) recently reported 
the discovery of OH from two HAeBe stars and concluded that the emission was a result of OH 
fluorescence. The rotational energy diagram of the transitions observed (in this case from the 
V = 1 state) indicate a single rotational temperature which leads the authors to conclude that it is 
unlikely to be the result of "prompt" emission, i.e., OH produced by the photodissociation of water. 
Their models of fluorescent emission from OH match the observed data from HAeBe stars quite 
well. In their models, OH is excited to the v=l state through rovibrational transitions pumped 
by near infrared photons and electronic transitions pumped by UV photons. However, the Spitzer 
data for TW Hya indicate mainly a downward rotational cascade from the high J states of v = 
state (Najita et al. 2010). Therefore, in TW Hya the observed OH is more likely to be prompt 
emission. 

H2 rotational lines Pure rotational S(2) and S(l) line emission from H2 was reported by Najita 
et al. (2010) at 4.8 x 10~^ and 10~^ Lq respectively. Most of the H2 emission in our models comes 
from within radii of 20 - 30AU, with ~ 35% of the H2 S(2) emission and 10% of the H2 S(l) 
emission from the inner SML disk (see Fig. [7]). In the mid-disk region, H2 emission originates in 
the superheated dust later at Ay ~ 0.3 — 0.1 and where the gas temperatures are 100 — 500K. H2 
formation in the mid and outer disk (r > 4AU) is mainly on grains and it is destroyed by FUV 
photodissociation. We obtain H2 S(2) and S(l) line luminosities to be 3.4 x 10~^ and 1.2 x 1O~^L0 
respectively. These are in reasonable agreement with the Spitzer IRS measurements. Bitner et al. 
(2008), in their ground-based survey for H2 using TEXES, did not detect H2 8(1) emission and 
obtained an upper limit of 5.0 x IO^^Lq. This suggests that half of the S(l) emission observed by 
Spitzer IRS must originate from a more extended region outside the narrower slitwidth of TEXES, 
i.e. at r > 20AU. Although the H2 S(l) emitting region is quite extended (Fig. [7|), half the H2 S(l) 
emission m the SML disk in fact originates at r < 14AU. The S(l) flux within 20 AU is 8 x lO'^L© 
and the radial extent of emission is in moderate agreement with the TEXES and Spitzer data. The 
line luminosity in the H2 S(0) transition is predicted from the SML disk to be 3 x 1O~^L0. 

[OI]6300A and 5577A emission We show below that the [OI]6300A line from TW Hya origi- 
nates from the disk, mainly from the photodissociation of OH (also see Storzer & Hollenbach 1998, 
Acke et al. 2005) with a small contribution from thermal emission. The line luminosities for the 
[OI]6300A and [OI]5577A hues from TW Hya are IO'^Lq and 1.4 x IO'^Lq respectively (Alencar 
& Batalha 2002, Pascucci et al. 2011, S.Edwards, private communication). We first show that the 
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observed [OI]6300A/[OI]5577A line ratio of '--^ 7 makes it unlikely that the origin is thermal. 

The oxygen atom has a triplet ground state (^P) and singlet D and S states which are at 
Eu/k ~ 22,850K and 48,660K respectively. The 01 6300A transition from the to the ground 
state and the 01 5577A transition from the state to the state therefore require high gas 
temperatures for collisional excitation to these high energy levels and subsequent radiative de- 
excitation. If the gas is collision dominated (in LTE), then the ratio of the two lines can be easily 
obtained from the ratio of UuA^E, where is the upper level population, A is the Einstein A- 
value and /S.E \s the energy of the photon emitted. In LTE, n{^S) ~ {l/b)n{''D)e^'^l^"^ , where 
(1/5) is the ratio of the statistical weights and /\E/k ~ 26, OOOK is the energy difference between 
the two levels. Therefore 

-^6300 ^ n(^Z:>)^6300A^6300 ^ 5 4 ^ iQ-'i^^mQO/T (^^^ 
-^^5577 'T-(^'S')^5577A-E'5577 

where ^6300 ~ 6 x 10~^ s~^ and ^5577 ~ 1.3 s~^. To obtain the observed line ratio of 7, gas 
temperatures therefore need to be ~ 3600K, which is possible with X-ray heating of the disk 
surface. However, the gas densities have to be high enough to produce LTE conditions. Collision 
partners can either be electrons or neutral hydrogen atoms and typical collision rate coefficients 
are 7e ~ 5 x 10^^(T/10^)'^'^ cm^ s~^ for electrons for both the transitions. The collisional rate for 
0(^D) with H atoms is 7// ~ 10"^^ cm^ s"^ (Krems, Jamieson &; Dalgarno 2006). "jh for 0(^S) is 
not known. However, the rates for the relaxation of the 0(^S) state are expected to be smaller than 
those for the relaxation of the state, because the spin-orbit couplings in electronically excited 
atoms are generally weaker (R. Krems, priv. comm.). We conservatively assume 7// ^ 10~^^ cm'^ 
s~^ for 0(^S) as well. Therefore, electrons dominate collisions if Xg > 4 x 10~^. For levels to achieve 
LTE, electron densities ng > A/^e ~ 10*^ cm~^ are needed for the 5577A line (rig > 10^ for the 
6300A line). Collisions with H atoms require neutral densities about 2500 times these values. Even 
for gas with high electron fractions Xg ~ 0.1, gas densities therefore need to be high, nu ^ lO^cm"^. 
These conditions are likely to be attained only in the inner disk, where gas at the surface may be 
dense, be partly ionized with sufficient neutral oxygen, and at T > 3600K. 

For gas at lower densities a balance between radiative de-excitation and collisional excitation 
and de-excitation of the different levels determine the level populations. We calculate the line 
ratios as a function of electron density and gas temperature and Figure [9] shows the contours for 
the observed value of 7 and other values about this number. From the contour plot, one can 
see that even for relatively high gas temperatures of ~ 8000K the required electron densities are 
~ lO^cm"^ implying gas densities of ~ lO^cm"^ for an electron fraction of Xg ~ 0.01, which is 
typical. Such conditions are unlikely, except perhaps in the innermost regions, r < 0.1 AU. We also 
note that soft X-rays that can heat the gas to such high temperatures typically have very small 
absorption columns ^ lO^^cm"'^, and very little mass is expected in gas at these high temperatures 
and densities. Such low masses may not be able to reproduce the observed high optical 01 line 
luminosities. 

For typical densities in the hot disk surface, n ~ 10^~'''cm~^, and electron fractions Xg ~ 0.01 — 
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10~^ the electron densities are rig < lO^cm"^, and the observed line ratio requires T > 30000K. We 
consider this extremely unlikely as it is very difficult to heat disk gas to these high temperatures 
and moreover retain a significant amount of neutral oxygen. The recent calculations of Ercolano 
& Owen (2010) also show that the 5577A line is appreciable relative to the 6300A line only when 
gas densities are high, as in their primordial disk models where the 01 emission is from the inner 
r < lAU region. In their models, the 01 emission is thermal in origin and they quote i630o/-f'5577 
~ 13 — 50 for the primordial disk models and ^ 300 — 700 for the transition disk models, supported 
by the above analysis. Only one of their many models, a primordial disk with Lx ~ 2 x 10^'^erg 
s~^ has a low line ratio ^ 13, where the emission presumably originates from gas with high electron 
density. 

In the SML disk, thermal emission contributes to about 20% of the observed [OI]6300A line 
and 10% of the 5577A line. In the inner disk, the luminosities are 1.2 x 10 ^L© and 10 ^L© for the 
thermal component of the 6300A and 5577A lines respectively, and the emission comes from regions 
where He ~ lO^-^cm'^ T ~ 4000 - 8500K at the disk surface. In the SMH disk, the calculated 
thermal [OI]6300A line luminosity is higher, 2 x 10~^Lq (TableE]). There is also thermal [OI]6300A 
emission from the mid and outer disk surface that is heated by X-rays to T > 5000K, the luminosity 
here is ~ 1.1 x IO^^Lq. There is negligible thermal [01] 5577A emission from the mid and outer 
disk regions. 

We propose here that most of the 01 forbidden line emission from TW Hya is the result of 
photodissociation of OH in the disk. The OH photodissociation layer in the SML disk is almost 
spatially co-existent with the H2O photodissociation layer, and this generally is located near the 
H/H2 transition in the disk. OH photodissociates about 50% of the time to produce an oxygen atom 
in the state, and occasionally (5% of the time) in the higher state (e.g., Festou & Feldman 
1981, van Dishoeck & Dalgarno 1983). Subsequent electronic de-excitation to the ground state 
produces the oxygen doublet at 6300A and 6363A (^D— ^P) and the ^S— transition produces the 
5577A emission. The atoms in the O^D state produce the 6300A-|-6363A doublet in the branching 
ratio of 3 (Storey & Zeippen 2000). The ratio of these line strengths is the same as the branching 
ratio because both lines originate from the same upper state. These calculations are supported by 
comet observations of the photodissociation products due to solar Lyman a (Morgenthaler et al. 
2007). We again count the number of photodissociations of OH in the model disk as we did for 
water, and multiply this by the energy of the transition and the theoretically calculated branching 
ratio. 

Our inclusion of the 0(^D) -|- H2 — > OH reaction leads to a lower effective photodissociation 
rate of OH as discussed earlier. If this reaction were not considered, then prompt emission produces 
a 6300A/5577A line ratio ~ 7. However, in regions where the densities are high, such as in the 
inner disk, the line ratio is expected to be lower as the reactive nature of the O^D atom will result 
m a suppression of 6300A emission. Although the state of oxygen is also more electrophilic than 
the triplet ground state, the unpaired electrons likely make the 0(^S)-|-H2 reaction less favorable. 
Further, the radiative decay of O^S is faster (^5577 > 200^630o) requiring even higher densities for 
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the reaction with H2 to be competitive. The 5577A hne is therefore not expected to be lowered 
significantly. In dense regions of disks, the ratio of 6300A/5577A is therefore expected to be lower 
than ~ 7. Nqh photodissociations per second of OH will produce ~ 0.5(1 — /o)-/V'ojy 6300A photons, 
where (1 — /o) is the reduction factor due to 0(^D) + H2 ^ OH which suppresses the radiative 
transition leading to the 6300A line. We therefore calculate that the [OI]6300A and [OI]5577A 
line luminosities due to OH photodissociation are 1.1 X IO^^Lq and 9.3 x IQ-'^Loin the inner disk 
respectively, and these lines originate in the same region as the OH cascade, at r < O.IAU, where 
n ~ 10^^~^^cm~^ and T ~ 500 — lOOOK. OH photodissociation in the mid-disk regions produces 
line luminosities 5 x IO-'^Lq (6300A) and 7 x 10-^Lq(5577A). This emission is produced from gas 
with n ~ 10^"^°cm~3 and T ~ 200 - 600K out to radii r ~ 40AU. Therefore, most of the 01 
luminosity from TW Hya comes from photodissociation of OH. Ercolano & Owen (2010) recently 
propose that the 01 forbidden line emission seen in disks may arise from a X-ray driven wind. Our 
picture differs from theirs in two ways - in our model the 01 emission arises through a non-thermal 
process, and the emission arises from the disk atmosphere rather than in a wind. This perspective 
is consistent with the 01 6300/5577A line ratio and the lack of a blueshift for the emission. 

The total luminosities from the SML disk (thermal -|-photodissociation) are therefore 8.4 x 
1O~^L0 for the [OI]6300A line and 1.6 x lO'^Lgfor the [OI]5577A line. Our 6300/5577A line ratio 
is ~ 6. The SMH disk which was a better match for the OH emision gives slightly higher [01] 
emission, 1.5 x lO^^L© for the 6300A hne and 1.92 x IO'^^Lq for the 5577A line for a ratio of 7.8 
(also see Table [6|). 

The intrinsic FWHM linewidth of the observed 01 lines are of order 10 km s~^ (Pascucci et 
al. 2011). This corresponds to a velocity dispersion of about 5 km s~^. In the SML disk, the 01 
lines arise mostly in the mid-disk regions where the projected Keplerian speeds are low, < 2 km 
s~^. In addition, the gas is relatively cool and turbulent speeds are expected to be low (Hughes et 
al. 2010). Hence, the contribution to velocity dispersion from thermal and turbulent velocity fields 
is expected to also be low, < 1 km s~^. However, the photodissociation of OH leads to a recoil of 
the resultant 0(^D) or 0(^S) atoms at speeds of order 2 km s~^, and the [01] photons are emitted 
prior to collisional relaxation (van Dishoeck & Dalgarno 1984). We estimate the linewidth of the 01 
6300A line by adding the various components and the Keplerian, thermal and recoil contributions 
as appropriate. We obtain a linewidth of ~ 6km s~^, narrower that the observed value. In the 
SMH disk, the calculated linewidth is more in agreement with observations, ~ 11km s~^, because 
of the increased prompt emission close to the star, although the total OI6300A line luminosity is 
overproduced by a factor of ~ 1.5 (Table [6|). 

We end with a caveat that although the above production mechanism via OH photodissociation 
explains the OI6300A and OI5577A line strengths, ratio, and lack of an observed blue-shift in TW 
Hya, it may not be a more general explanation valid for all disks in which optical forbidden line 
emission has been detected, especially if these lines are blue-shifted and originate in a flow. For 
disks with different surface density distributions and FUV and X-ray luminosities, there may be a 
more dominant thermal contribution to [01] emission. Further, our disk models are static and do 
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not include any emission that might arise in an extended wind (e.g., Ercolano and Owen 2010). 
Although we expect that our disk density and temperature structure is reasonably accurate within 
our grid, we do not account for any emission that might arise from heights above our vertical 
extent of z = r. Densities in these regions are typically low, and usually not a significant source 
of emission. For the optical lines (01 and SII, discussed below) however, the upper states are at 
~ 20, 000 — 40, OOOK and the emission is very sensitive to temperature, and the contribution from 
the wind in this case might be significant. 

SII optical lines Pascucci et al. (2011) re-analyze the optical data of Alencar & Batalha (2002) 
and report the detection of the [SII]4069A line (~ 3 x 10~^Lq) from TW Hya. They do not detect 
other lines of SII and place an upper limit of 1.7 x 1O~^L0 on the 673lA line. We find that in our 
SML disk, most of the SII emission arises from the inner disk in hot, dense gas that also results in 
the thermal 01 emission calculated above. Our calculated SII line luminosities are 3.9 x 1O~^L0 
for the 4069A lines, and 8.6 x lO''^, 5.6 x 10"^ and 1.1 x W'^Lq for the undetected lines at 4076A, 
6716A and 673lA respectively. Our [SII]4069A line is slightly stronger than observed, however, 
sulfur may be somewhat depleted compared to our assumed abundance(2.8 x 10~^). 

Nell and Nelll [Nell] emission from TW Hya is strong and has been detected by Ratzka et 
al. (2007), Herczeg et al. (2007), Pascucci & Sterzik (2009) and by Najita et al. (2010). There 
is some discrepancy between the measured values, with more flux observed by the Spitzer study. 
Najita et al. (2010) conclude that [Nell] emission arising from a more extended region (r > lOAU) 
explains their higher flux at all three observational epochs, which the narrow slits of the ground- 
based observations exclude (Pascucci & Sterzik 2009). A revised calibration (Pascucci et al. 2011) 
recovers the Spitzer flux, and suggests that most of the [Nell] emission comes from within ~ lOAU. 
Earlier work by Herczeg et al. (2007) measured a line luminosity of 4.4 x lO~^L0with their ground- 
based study, similar to the Spitzer data. There is therefore some indication that the [Nell] line flux 
(or the underlying continuum) might be variable (Pascucci et al. 2011). Here we use the Spitzer 
IRS measurements of Najita et al. (2010) for the [Nell] and [Nelll] emission, 3.9 - 4.7 x lO'^L© 
and 2 x IO'^Lq. 

[Nell] emission from the disk models is found to be quite insensitive to the surface density 
distribution. This is because in X-ray heated regions, the [Nell] emission mainly depends only 
on the intercepted X-ray photon flux and the temperature of the emitting gas (HG09). A simple 
estimate of the expected [Nell] line luminosity can be made using the analytical expression derived 



where (px ~ lO^^s"^ is the IkeV X-ray photon luminosity. In the above, we also assume that the 
disk intercepts a fraction / ~ 0.7 of the stellar X-ray flux (our grid in the numerical work typically 
extends to Zmax = r). A typical gas temperature in the [Nell] emitting layer is T 2000K. The 
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expected [Nell] luminosity for the TW Hya disk is therefore, 4 x 10 ^Lq. 

From our full, detailed numerical calculations, we obtain total [Nell] and [Nelll] line luminosi- 
ties for the SML disk as 3 x 10~^Lq and 3 x 10~'^Lq respectively. Approximately 25% of the [Nell] 
comes from gas in the inner disk (see Fig. [7]). The mid-disk region (3.5AU < r < 20AU) contributes 
the remaining 75% to the [Nell] luminosity in the SML disk. The [Nell] arises from the X-ray heated 
and partially ionized surface of the disk where gas temperatures range from 1000 — 4000K. [Nell] 
emission from the SML disk is restricted radially to r < 10 AU (Fig. [7|), and peaks strongly at the 
vertically extended 4 AU rim. The X-ray produced [Nell] is somewhat lower than the observed line 
flux, which could be explained either by the presence of an EUV component and/or by variability 
of the stellar high-energy photon flux. We note that the accretion rate of TW Hya is observed to 
be variable on timescales similar to that for the [Nell] line flux (Eisner et al. 2010), ~ a few years. 
As the X-ray and any [Nell] producing EUV radiation fields may be a result of accretion-induced 
stellar activity (Robrade &: Schmitt 2006, HG09), this may naturally explain the variability of the 
[Nell] line. 

Pascucci & Sterzik (2009) resolve [Nell] in their spectra, and the observed linewidths and 
blue-shifted emission are consistent with a photoevaporative flow from the star. The line profile is 
very well reproduced by EUV {hiy >13.6eV) photon-driven photoevaporation models (Alexander 
2006). However, we do not include EUV in our source spectrum for TW Hya and can recover 
most of the observed flux from X-ray heated gas alone. The remaining ~25% of the observed 
Nell flux may be due to the EUV layer. We can thus estimate the EUV photon luminosity as 
being ^ 3 x lO'^'^ s~^. The [NeII]/[NeIII] line ratio may provide some measure of the origin of the 
emitting gas (the completely ionized EUV-heated layer or the partially ionized X-ray-heated layer 
below the EUV-heated layer; 10,000K in the EUV case and ~ 1,000 - 3,000K for X-rays), but is 
not definitive because of the uncertain nature of the EUV spectrum (HG09). However, the most 
natural explanation is an origin in the X-ray heated gas, where rapid charge exchange of Ne"*"^ 
with H atoms quenches the [NeHI] line and produces the observed [Nen]/[Neni] ratios. (For a 
more detailed discussion, see HG09.) The observed blue-shift may also arise for a X-ray induced 
photoevaporative flow (Gorti et al. 2009, Owen et al. 2010), and the lack of need for significant 
EUV to reproduce the [NeH] suggests that any photoevaporative flows may originate in primarily 
neutral. X-ray heated gas. 

The far-infrared and sub-millimeter emission lines discussed hereon originate only in the mid 
and outer-disk regions and are the same for the two disk models, SML and SMH. The line lumi- 
nosities are listed in Table [6l 

[OI] and [CII] fine structure lines We obtain a [OI]63/um line emission of 6 x 1O~^L0, a factor 
of ~ 2 higher than the Herschel PACS measurement of Thi et al. (2010). Our calculated [OI]145^m 
line luminosity is 3 x 1O~'''L0, and consistent with the upper limit of 4 x 1O~'''L0 obtained by Thi 
et al. (2010). Emission is extended and produced by ~ 50 — lOOK gas at r ~ 30 — 120AU. 
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Our model overproduces [OI]63//m emission, and we suggest that inclusion of water freezing 
on grains might alleviate this problem. Freezing of water ice on grains colder than lOOK removes 
gas phase elemental oxygen from the outer disk where dust grains are cold. Photodesorption would 
release some of this ice and at Ay < 3, all water ice is photodesorbed from grains (Hollenbach et 
al. 2009). Including these processes reduces the [OI]63/im emission from the disk compared to a 
disk with no ice formation, to 4 x 1O~^L0, in better agreement with the PACS observations. The 
[OI[145/Lim line is reduced from 3 x 10~^Lq to 2 x 1O~^L0. A more detailed treatment of freeze-out 
of molecules on dust grains and desorption processes will be a subject of future work. Another 
possibility is that the PAH abundance in the outer disk is even lower than we assume, with a 
depletion factor of nearly ~ 1000. This lowers the temperature in the outer disk sufficiently and 
decreases the [OI]63^m line luminosity to 4 x 1O~^L0. If freezing were also included, this would 
bring the 01 emission even closer to the observed value. However, we note that decreasing the PAH 
abundance everywhere in the disk alters the vertical structure of the mid-disk region (r < 20AU) 
sufficiently to deteriorate fits to observation of other species such as H2, OH and 01 forbidden lines. 
There may, in fact, be a radial gradient in the PAH abundance in the disk, with a higher degree of 
depletion in the outer disk, if PAHs were to enter the disk frozen out on dust grains as suggested 
by Geers et al. (2009). Another possibility is that the PAH abundance is low everywhere in the 
disk, but that other sources of heating such as mechanical viscous dissipation are effective at the 
surface in the disk(r < 30AU) (e.g., Glassgold et al. 2004). 

Thi et al. (2010) model [01] emission from TW Hya and conclude that the observed flux 
together with their ^^CO observations indicate a low gas disk mass by factors of 10 — 100 compared 
to that expected from the disk dust mass and assuming a gas/dust ratio of 100. However, we find 
that lowering the gas disk mass by 10 and retaining the dust mass only marginally improves model 
fits to the observed [01] emission, and worsens the other line emission fits considerably. Most of 
these discrepancies can be attributed to the increased grain absorption cross section per H in a disk 
with a low gas/dust ratio. H2 emission increases (FUV heating by PAHs becomes significant) by a 
factor of ~ 2 — 3, OH prompt emission is reduced due to the increased (relative to gas) absorption 
of Lyman a photons by dust, and there is excess emission from CO pure rotational lines. The 
HCO"''4-3 line requires high densities for emission, and this line luminosity is almost a factor of 10 
lower for the low mass disk than the observed value. Increasing the PAH abundance by 10 as in 
the preferred model of Thi et al. would raise gas temperatures even higher and further increase H2 
and CO emission. Although we do not concur with the Thi et al. (2010) results, we note that our 
calculated model line fluxes are in good agreement with results from other models they consider. 
Our model is most similar to their low PAH abundance models (their series 2 and 3) and if they were 
to base their conclusions on these alternate models, a disk mass close to ~ O.O6M0 is inferred from 
their results. For these series, they also seem to find that the 01 line luminosity is quite insensitive 
to disk mass in the range we consider. Their conclusion of a low gas disk mass appears to be largely 
based on the ratio of the ^^CO/^^CO line fluxes. In order to accurately interpret this measurement, 
a more accurate disk chemical model which includes ^'^CO as a species is necessary. Our model 
does not include ^^CO. We also note that Thi et al. assume a constant ^^CO to ^^CO ratio to 



- 26 - 



calculate the CO flux and do not include a proper chemical model of CO as well. Moreover, the 
^^CO/^^CO line ratio will be significantly affected by processes such as selective photodissociation 
(Lyons &: Young 2005) and freeze-out processes, which are not included in our disk models or those 
of Thi et al. (2010). It remains to be seen whether a low gas disk mass is needed to reproduce the 
observed CO isotope flux ratio. 

The calculated [CII]158//m line luminosity is 3 x 10~'^Lq, and is consistent with the upper 
limit of 4.5 x 1O~^L0 set by the Herschel observations. 

CO rotational lines The observed line luminosities of the CO 6-5, 3-2 and 2-1 lines are 5.0 x 
10^*^,2.5 X 10"^ and 7.9 x 1O~^L0 respectively, whereas we obtain 8 x 10~^ 2.4 x 10"^ and 
7.9 X IO^^Lq for these lines. Purely thermal line broadening in the outer disk provides better 
correspondence between data and model results. This is in fair agreement with the results of Qi et 
al. (2006) who assume only a small amount of turbulence (Ay ~ 0.3 kms~^) in their analysis and 
with the more recent results of Hughes et al. (2010) who find that the turbulent linewidths are 
low. A power-law surface density profile also matches the data, but has to be arbitrarily truncated 
at ~ 70AU (GH08). We consider it more realistic to assume that the surface density declines more 
smoothly with radius, as in a disk evolving viscously and subject to photoevaporation (also see 
Hughes et al. 2008, Isella at al. 2009). 

CO rotational line emission constrains the distribution of gas in the outer disk. We do not 
model disk masses higher than 0.06Mq, because we expect that appreciably higher disk masses will 
make the disk gravitationally unstable. Lower disk masses by a factor of 10 (Vldisk = 0.006Mq) 
do not fit the CO rotational line emission well. We find that although it is possible to match the 
observed CO 2-1 and 3-2 line emission by altering the surface density profile in the outer disk for 
the lower mass disk, the CO 6-5 line emission is higher by a factor of ~ 4 than the SML value, 
and hence discrepant with data. We expect freezing will not affect the emission line luminosities 
calculated. A simple estimate made by discounting contributions from spatial locations where 
the dust temperature is < 20K only slightly lowers the calculated line luminosities by ~ 15%. 
Photodesorption may also keep the CO molecules from freezing on dust grains in the surface 
regions of the outer disk. We note that we do not need severe radial depletion of CO to explain 
the observed line ratios (e.g., Qi et al. 2006), although some degree of depletion (less than a factor 
of 3) in the mid-disk region may help lower the CO 6-5 line luminosity from the model disk. 

HCO^ 4-3 line van Zadelhoff et al. (2001) measured the line intensity of the HCO"*" 4-3 transition 
from TW Hya using the JCMT, and the line luminosity is 7.7 x 10~^Lq. Our calculation yields 
4 X 1O~^L0, nearly a factor of 2 lower than observed. We could obtain a better match by changing 
the surface density distribution in the outer disk to increase the strength of the HC0"'"4 — 3 line, 
but this then also changes the CO rotational line intensities and makes that agreement poorer. 
Further, increasing the outer disk surface density also results in a high overall disk mass, bringing 
the disk closer to being gravitationally unstable. The HCO^ line mainly originates in relatively 
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warm (T ~ 20 — lOOK), and high density gas (n > 10^~^cm~^) at r < 50AU. The abundances 
of HCO"'' calculated by the model in the emitting layers of the disk are a factor of 20 — 70 lower 
than the typical dark cloud abundance of 5 x 10~^ adopted by van Zadelhoff et al. (2001) in their 
analysis, and thus consistent with their derived "depletion" factors of 10 — 100. 

We conclude this section with the caveat that our models may not be unique solutions to the 
gas surface density distribution in the disk. There are theoretical uncertainties that are inherent in 
chemical rates and the microphysics, as well as several important processes (such as freezing of ices) 
have been either simplified or ignored. The complicated and computationally intensive modeling 
procedure also makes a full exploration of parameter space impossible at this time. Nevertheless, 
we find it promising that theoretical models with reasonable parameters can reproduce observed 
emission ranging from the optical to the sub-millimeter wavelengths with an accuracy to within a 
factor < 2. 

4. The Origin of the Inner Gas Hole 

Modeling of gas emission lines and comparisons with observational data indicate that the 
surface density of gas is also depleted in the inner r < 4AU dust depleted region. The models are 
further consistent with a massive, optically thick outer disk beyond the 4AU rim. At least three 
different explanations have been proposed to explain the dust holes in transition disks, of which 
TW Hya is a leading example (e.g., Najita et al. 2007). The distribution of gas in the disk that we 
derive may help distinguish between these various scenarios. We examine the implications of each 
of the explanations for transition disks and discuss our model results in each context. 

Grain growth with no gas depletion in the dust hole: The first step towards planet forma- 
tion in disks via the core accretion scenario is the coagulation of primordial sub-micron sized dust 
grains into larger dust particles and planetesimals. Grain growth is expected to be more efficient 
in the inner disk as the densities there are higher (e.g., Lissauer et al. 2007). The population of 
small dust decreases due to grain growth, and as larger solids are not efficient at absorbing starlight 
because of their lower collective surface area, this process is manifested as a decrease in the dust 
continuum emission at shorter wavelengths where the warm small dust in the inner disk would 
emit. Grain growth does not remove disk gas, and no depletion in gas surface density within the 
inner regions is expected (unless gas giant planets have already formed, which is the next scenario). 
In this case, there would be continued accretion onto the star (as is observed for TW Hya). 

Our modeling of the gas disk around TW Hya when compared with observational data does 
not support the grain growth explanation alone as a mechanism to cause the dust hole. Based on 
the constraints set by gas line and continuum emission, we argue that gas must also be depleted 
in the inner disk region (r < 4AU), with a surface density decrease ~ 1 — 2 orders of magnitude 
compared to what one would obtain by extrapolating the surface density from the outer optically 
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thick, massive disk. Although grain growth may explain the decrease in dust opacity inside 4 AU, 
a mechanism that depletes gas surface density by ~ 10 — 100 is also necessary. 

Giant Planet Formation A widely favored explanation for the observed dust depletion in tran- 
sition disks is the creation of a gap by a giant planet, followed by the formation of a hole as the 
inner disk accretes onto the star. A gas giant planet exerts dynamical torques on the disk and 
clears material in its immediate vicinity to create a gap (e.g. Artymowicz & Lubow 1996, Bate et 
al. 2003, Rice et al. 2003). The disk interior to the planet is expected to rapidly accrete onto the 
central star. The outer disk mainly accretes onto the planet, with a small fraction of the mass flow 
(depending on the ratio of planet mass to disk mass) continuing to accrete onto the star (Lubow 
and D'Angelo 2006). Pressure maxima in the disk just outside the planet may act to filter dust 
particles (Rice et al. 2006) and further deplete the inner disk of dust. Our modeled surface den- 
sity distribution supports this explanation for the TW Hya disk, with the potential planet located 
between 2 — 4AU where the surface density may be quite low (with gas in accretion streams) and 
consistent with the presence of a gap. The accretion streams then feed a depleted inner accretion 
disk at < 2AU. 

We can use the surface density distribution determined by the gas emission models to set 
constraints on the mass of the planet orbiting TW Hya. As the planet must open a gap in the disk, 
or at least establish a large surface density contrast at its location, we can use the gap-opening 
criterion to obtain a minimum mass Mp of the planet (e.g., Lin &: Papaloizou 1986) 



where H is the scale height of the outer disk at the radius r ~ 4 AU. From our disk models, H/r ~ 
0.04 at 4AU, and therefore the planet has a minimum mass of ~ 0.03 — 0.3Mj for a ~ 10"'^ — 10~^ 
in the outer disk. Simultaneously, one can appeal to recent hydrodynamical models of accretion 
onto and past a massive planet (Lubow &: D'Angelo 2006), which find that a planet with at least 
this minimum mass must accrete matter at a rate ~ 10 times higher than the rate onto the star. 
Given the current age of the system (~10 Myrs), the planet is likely to have undergone accretion 
for at least 1 Myr. Even with the lower estimate of the stellar accretion rate ~ 4 x IO^-'^'^Mq yr~^, 
this implies a minimum planet mass Mp ~ 4Mj. An upper limit to the planet mass can also be 
estimated from the stellar accretion rate. The fact that gas in the outer disk can still flow through 
the gap induced by the planet suggests the planet is not massive enough to form a tidal barrier 
and repel disk matter. Depending on a, this limits the planet mass to ~ 7Mj{a = 4 x 10~^) or 
~ 2Mj{a = 4 X 10^^) (see D'Angelo et al. 2010). We conclude that a plausible planet mass is 
likely to be ~ 4 - 7Mj. 

The mass of the planet inferred above is high enough that dust particles with sizes larger than 
lyLim may be trapped in the pressure maximum in the disk at the outer edge of the gap (Rice et 
al. 2006). Dust particles with sizes of ~ 2/xm have been modeled to exist in the inner disk (C02), 




(3) 
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consistent with this scenario. Since the inferred dust distribution in the outer disk has a maximum 
size of a^nax = 1cm, and only sub-micron sized grains cross the gap, dust filtering would result in a 
decrease in dust mass inside the gap. Therefore an increase in the gas/dust mass ratio across the 
gap by a factor of ~ 100 (since the dust mass is proportional to ^/{amax)) is expected, which is 
consistent with our SMH model disk. 

A potential issue with the planet hypothesis is the migration rate. For this case, we need to 
consider Type II migration as the planet must be massive enough to create a gap. Viscous diffusion 
of the disk causes the planet to drift towards the star efficiently, but only as long as the local disk 
mass (7rr^5](r), r 4AU) is greater than the planet mass. For the mid and outer disk surface 
density distribution, the local disk mass is ~ IMj. Thus the expected planet of ~ 4Mj is stable 
against migration. 

Photoevaporation 

The strong FUV and X-ray flux from TW Hya should drive massive photoevaporative flows from 
the disk (Ercolano et al. 2008, Gorti k Hollenbach 2009, Gorti et al. 2009, Owen et al. 2010). 
Photoevaporating flows have been observed by resolved ground-based [NeII]12.8//m emission from 
TW Hya (Pascucci &: Sterzik 2009). The observed blue-shifted line proflles are well-reproduced by 
theory (Alexander 2006, Owen & Ercolano 2010). 

However, whether the inner hole in TW Hya is carved by photoevaporation remains a question. 
The location of the radius of the hole ~ 4 AU is consistent with theory as it is larger than the 
so-called critical radius, ~ lAU where the hole flrst appears in photoevaporative theory before 
the disk is eroded radially outward by direct rim illumination. However, photoevaporation theory 
traditionally predicts that a "clean" hole, completely devoid of dust and gas, must rapidly form 
on viscous timescales (e.g., Alexander et al. 2006). The viscous clearing timescale in the inner 
disk is tyis ~ r'^/u, or t^is ~ y^(GM*r)/(aCg). For the TW Hya disk, t^s ~ 10^ years at 4AU 
for a ~ 0.01. Therefore, the surface density drop of ~ 100 at 4 AU implies an e-folding factor of 
5, or ~ 5 X 10^ years after gap opening by photoevaporation. For a higher surface density in the 
inner disk and a drop of ~ 10 at 4 AU, gap opening must be even more recent. Although there 
must be ongoing photoevaporation from the disk as is observed, our results appear to refute current 
photoevaporation theories for the creation of the inner hole in TW Hya, unless the gap creation 
epoch is extremely recent, which seems unlikely given the ~ 10^ year age of TW Hya. 

We note that one possibility where one could avoid a special observing epoch and where 
photoevaporation could work as a gap creation mechanism would be if the thermal photoevaporative 
flow that emanates from the rim and disk surface is partially re-captured by stellar gravity to form 
the inner disk and eventually accrete onto the star. In our models, the photoevaporative mass loss 
rate is ~ 4 x 10~^MQyr~^. For a stellar accretion rate of ~ 4 x 1O~^^M0 yr~^, 10% of the flow needs 
to be recaptured by the star. We use simple timescale arguments to show that this is a feasible 
process. For recapture, the flow timescale (tj) to the sonic radius where gas escapes needs to 
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be comparable to the viscous timescale tyis of the photoevaporating gas. As tj ~ Ts/cg, and using 
'i^s = i^gf^ and the expression for t^s in the preceding paragraph, we obtain tf/tms ~ a{rg/r)^/'^ . 
For FUV/X-ray driven photoevaporation where the flow is launched subsonically (GH09, GDH09) 
from the rim at r = 4AU, and a typical gas temperature of ^ 2000K, rg = 31AU. For 10% of the 
flow is to be recaptured due to viscosity, we need tf/t^is 0.1 or a value of the viscous parameter 
a = 0.1. If the accretion rate is ^ 10"^Moyr-\ the mass captured needs to be ~ 25% for the 
calculated photoevaporation rate. Therefore, for this re-capture process to be possible, a needs to 
be relatively high in the inner disk, ~ 0.1—0.25. However, this process is a complicated phenomenon 
that requires full 2D hydrodynamical calculations to resolve, including turbulent viscosity in the 
flow that can redistribute the angular momentum. This problem will be the subject of future study. 

We also note that even if the inner opacity hole is not created by photoevaporation, the 
process must still act to disperse the disk. The observed [Nell] line profiles indicate this is an 
ongoing process in the TW Hya disk. We can estimate the lifetime of the remaining disk around 
TW Hya using FUV/X-ray photoevaporation rates from theory (GH09) and the results from the 
SML disk. From our calculated disk photoevaporation rate 4 x 1O~^M0 yr~^, we derive a disk 
lifetime of M^isk/Mpe as ~ 10^ years. Note that with the reduced disk mass of Thi et al., the 
expected disk lifetime is < lO^years, which places us in an unlikely viewing epoch and therefore 
seems implausible. We calculated a time evolution model of the disk, with the SML disk as an 
initial configuration and assuming a viscosity parameter a = 0.005 in the outer disk using the 
photoevaporation models described in Gorti et al. (2009). When the effects of viscous evolution 
are included, the remaining lifetime of the disk is calculated to be 5 x 10^ years. We note that a 
higher assumed value of a would lead to even shorter disk lifetimes. 

Most of the photoevaporation occurs at the 4AU rim, although the mass loss rate derived is 
lower compared to that from holed disks with directly illuminated rims (e.g., Gorti et al. 2009, 
Owen et al. 2010) because of the partial shielding of the rim in this case by gas in the inner disk. 
The photoevaporation rate is higher than recent measurements of the stellar accretion onto TW 
Hya, which range from 4 X 10~^° - 3 x lO'^M© yr~^ (Muzerolle et al. 2000, Alencar & Batalha 
2002, Eisner et al. 2010). This is consistent with the requirement for photoevaporation to create a 
gap in the disk, when the accretion rate at ~ 1 — 4AU must drop below the photoevaporative mass 
loss rate. 

If photoevaporation has caused the inner hole and not a giant planet, then the amount of gas 
we infer to be present in the inner disk and the observed accretion onto the star imply that either 
the gap has been very recently created (< 10^ years ago), or, that the gap was created earlier and 
that the observed accretion is caused by viscosity in the fiow causing significant re-capture of the 
photoevaporating gas and accretion onto the star. If a giant planet has caused the inner hole, then 
the accretion rate in the outer disk is expected to be at least a factor of 10 higher than M*, and 
in this case the derived photoevaporative mass loss rate is lower than the outer accretion rate, 
which is consistent with the giant planet hypothesis and, in this case, not consistent with the pure 
photoevaporation scenario for hole formation. 
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5. Constraints on the Disk Mass from Line Emission Modeling 

In this paper, we model the gas disk of TW Hya and conclude that the best model which 
reproduces the observed gas emission line luminosities is one where there is an optically thick 
disk (O.O6M0) extending from ~3.5AU outwards, with a depleted inner region that contains ~ 
10~^ — 10~^Mq of gas extending to the star. The modeling suggests a lower gas mass in the inner 
disk by ~ 10 — 100 than obtained by extrapolating the surface density of the outer disk radially 
inward. Most of the disk mass (O.OGMq) is contained within lOOAU, and the surface density in 
the outer disk declines sharply beyond this radius. The gas/dust mass ratio is consistent with the 
canonical ratio of 100 in the outer optically thick disk, and is enhanced by factors of ~ 5 — 50 in the 
inner disk where both gas and dust are severely depleted compared with the extrapolated surface 
density of the outer disk. Although the outer disk gas/dust ratio is the same as the interstellar 
ratio, the cross sectional area of dust per H nucleus is 100 times less than interstellar dust, 
because the disk dust size distribution extends to much larger sizes, Umax ^ 1mm, than interstellar 
dust where amax ~ 0.2//m. Therefore there are far fewer small dust particles (which provide most 
of the cross sectional area) per H nucleus, since much of the dust mass is tied up in the largest size 
(~ 1 mm) particles. 

Inner Disk The observed CO rovibrational emission and limits on H2O line emission between 
10-19 /im constrain the inner disk mass, and are best explained by a gas mass ~ 10"^ Mq within 
r < 3.5AU. Lower disk masses (1O~^M0) are very discrepant with observed high J CO rovibrational 
emission and can be ruled out. A full, undepleted disk with mass ~ 1O~^M0 is rejected based on the 
expected continuum emission due to gas opacity and the lack of observed mid-IR water emission. A 
disk mass of 10~^Mq is, however, also in reasonable agreement with data. There is some uncertainty 
in the inner disk dust mass as inferred from observations (~ 2 x 10~^Mq,C02; ~ 10~^'^Mq, Eisner 
et al. 2006), but the conclusion that the hole is depleted of gas and dust compared to the outer 
disk appears quite robust. 

The amount of mass present in the inner disk is fairly consistent with the observed accretion 
rate. Assuming steady state accretion rate onto the star M^^ = Svrz/E, where the viscosity u = ac^ /Q, 
a being the viscosity parameter, Cg the sound speed of gas and the angular velocity (e.g., Pringle 
1981), we can estimate a from the observed and from S and Cs derived from our models. For 
an accretion rate M* = 4 x 1O~"^''M0 yr~^, and at r = lAU, we then obtain a ~ 0.04 for the SML 
disk and ~ 0.004 for the SMH disk. (If were ~ IO^^Mq yr-\ then a = 0.1 for the SML disk 
and a = 0.01 for the SMH disk.) At the inner edge, close to the magnetospheric truncation radius 
of 0.06AU, a is higher, ~ 0.08 — 0.008. Our derived values for a for the preferred SML disk are 
somewhat higher than the 0.01 — 0.005 value typically used for protostellar disks, but consistent 
with the higher gas temperatures and ionization levels expected in the inner disk ( King et al. 2007, 
Turner et al. 2010). For the SMH inner disk mass of 1O~^M0, the inferred value of a is more in 
line with values typically used. 
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Outer Disk HCO"*" and CO rotational lines place upper limits on the amount of gas that is 
present in the outer disk, and are consistent with a disk gas mass ~ 0.06Mq. All the other 
observed line emission originates mainly in the disk at radii r < 30AU, which contains about 30% 
of the total disk mass. We find that decreasing the gas mass by a factor of 10 lowers the HCO"*" 
luminosity by a factor of 5 due to its origin in dense gas. Since our model already underestimates 
this luminosity by a factor of 2, this would produce a major discrepancy in the model luminosity 
versus the observed HCO"'' luminosity. In the adopted surface distribution with a tapering outer 
edge, CO pure rotational lines are also quite sensitive to the mass in the outer disk. The lower 
mass disk results in a high CO 6-5/2-1 line ratio of 16, compared with the observed ratio of 5.8 
(see Table [1]). Increasing the gas mass by a factor of 10 brings the disk mass to nearly the stellar 
mass, and such a disk would be highly gravitationally stable (Pringle 1981). We therefore consider 
a more massive disk than O.O6M0 very unlikely. In the outer disk, gas line emission may therefore 
constrain disk gas mass to within a factor of < 10, even though most of the gas line emission 
originates in the surface layers which do not trace the bulk of the disk mass. 

We can also set limits on the disk mass from arguments based on the presence of a planet and 
from expected mass loss rates due to photoevaporation. If a planet were present and accreting at 
~ 10 times the stellar accretion rate, this would imply a relatively short disk survival time for a 
O.OO6M0 disk. Even for a low accretion rate onto the star of M* = 4 x 1O~^'^M0 yr~^, the disk 
survival time would be ~ 6 x lO~'^M0/(lOAf*) ~ IMyr, which is short compared to the age of the 
system. This situation is made even more dire if the disk mass were lower or if the stellar accretion 
rate is at the higher range of measurements. Further, if the disk is losing mass at the estimated 
photoevaporation rate of 4 x 1O~^M0 yr~^, a O.OO6M0 disk would dissipate in < IMyr, which seems 
unlikely given the age of TW Hya. In view of all the above, our preferred gas disk mass is O.O6M0. 

6. Conclusions 

We infer the gas surface density distribution of the transition disk around TW Hya by modeling 
gas line emission from the disk, including H2, HCO"*" and CO rotational lines, CO rovibrational 
lines, OH line emission, [Nell], [Nelll], [01] fine structure emission and optical 01 and SII forbidden 
line emission. Our standard best fit model (SML) disk reproduces observed line emission well to 
within a factor of ~ 2 for most lines, except OH emission which is discrepant by a factor of ~ 4. 
We consider prompt emission following photodissociation of water as mechanism to produce OH. 
If we can increase our gas phase water production at the disk surface by a factor of a few, we can 
bring all emission lines in the model to within a factor of 2 of the observations. Vertical advection 
of icy grains in the mid disk region to surface layers where the ice mantles are photodesorbed may 
be the most likely source of the extra water vapor. We further speculate that photodissociation 
of OH itself may result in re-formation to a rotationally excited state in very dense regions, which 
might lead to the observed OH emission. The alternate model of a higher inner disk mass (SMH) 
provides good matches to the observed emission as well. The SMH disk is better at reproducing the 
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OH emission (prompt) and the [OI]6300A linewidth, although the CO rovibrational hne emission 
is a poorer fit and some emission lines are over-produced. Our main results on gas line emission 
are as follows. 

1. The model disk that best matches the observational data indicates a significant decrease in gas 
surface density in the inner opacity hole that causes a dip in the mid-infrared continuum. Our 
best models have a gas mass of 10~^ — 10~^Mq at 0.06AU < r < 3.5AU, and O.OGMQbetween 
3.5AU < r < 200AU. We exclude a full undepleted disk with 1O~'^M0 and inner disk masses 
of 1O~^M0 and below. The gas/dust mass ratio in the inner disk is inferred to be greater by 
a factor of ~ 5 — 50 than the canonical value of 100 which applies in the outer disk. 

2. CO rovibrational emission and upper limits on H2O rotational lines provide good constraints 
on the mass in the inner disk. The observed CO rovibrational emission is modeled to arise 
from hot (T ~ 500 — lOOOK) gas inside ~ 0.5AU, from disk surface regions where hydrogen 
is predominantly atomic. 

3. The mass of the outer disk is constrained to lie within a factor of a few of 0.06 Mq. Sig- 
nificantly higher masses lead to gravitational instability. Significantly lower masses lead to 
very poor fits to the CO rotational spectrum and to the HCO"*^ J= 4-3 line luminosity. In 
addition, given the relatively rapid accretion and photoevaporation of the outer disk, lower 
disk masses give unreasonably short outer disk lifetimes (< IMyr). 

4. We propose that OH mid-infrared emission originates from prompt emission following pho- 
todissociation of water in the disk very close to the star (r < 0.1 AU) and from the mid-disk 
regions (r ~ 4 — 40AU). Our model underestimates the production of OH by a factor of 4, 
and we speculate that vertical mixing may bring more water to the surface to be photodis- 
sociated. Grain formation of water ice followed by photodesorption could also increase the 
water vapor at the surface and thereby also increase the OH emission. However, enhanced 
water production might also lead to enhanced [01] 6300 A forbidden line emission by a factor 
of 3 above the observed value. We suggest another promising mechanism for producing 
OH emission in dense regions where the photodissociation of ground state OH leads to an 
excited O^D atom that reacts with H2 to form rotationally excited OH that cascades giving 
rise to the observed MIR lines. 

5. We argue that the forbidden lines of [01] at 6300A and 5577A arise due to the photodissocia- 
tion of OH mainly from the mid-disk regions (r ~ 4 — 20AU). Such an origin can successfully 
explain the OI6300A/5577A ratio of ~ 7 observed. About ~ 20% of the [OI]6300A emission 
is thermal in nature and has equal contributions from the inner disk (r < 4AU) and the outer 
disk (r ~ 4 - lOAU). 

6. Nell emission may not require the presence of a significant EUV flux from TW Hya and can 
be explained by X-ray heating of gas. If EUV is present, we estimate a photon luminosity of 
< 3 X lO^'^s^^ in order to explain the slight deficit in X-ray produced Nell line fluxes. 
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7. H2 pure rotational line emission is extended and is mainly from radii r 4 — 30AU. [OI]63/xm 
emission, CO pure rotational lines and HC0"''4— 3 line trace the outer disk at r 30 — 120AU. 

We conclude that grain growth alone cannot be responsible for the dust hole in the TW Hya disk 
because gas line emission modeling indicates a factor of ~ 10 — 100 depletion in the gas surface 
density as well. The transition disk morphology of TW Hya may be explained by the presence of a 
~ 4 — 7Mj planet at r ~ 2 — 3AU that causes the change in surface density distribution near this 
region. The disk is presently photoevaporating with an estimated mass loss rate of 4x 10~^MQyr~^. 
If there were no planet, then photoevaporation alone may be in the process of developing a hole. 
The calculated photoevaporation rate is greater than the stellar accretion rate suggesting that 
photoevaporation could have created the gap. However ongoing stellar accretion at current rates 
poses problems for present photoevaporation theories. Regardless of the hole-creation mechanism, 
the remaining disk lifetime due to photoevaporation and viscous evolution is estimated to be ~ 5 
million years. 
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Table 1: Observations of gas emission lines 



Line 


A 


Line Luminosity (1q) 


Reference 


Obs. Facility 


CO 6-5 


434/im 


4.6 X 10"^ 


Qi et al. 2006 


SMA 


CO 3-2 


867/im 


2.4 X 10"^ 


Qi et al. 2006 


SMA 


CO 2-1 


1.3 mm 


7.9 X 10"^ 


Qi et al. 2006 


SMA 


HCO+4-3 


841/im 


7.7 X 10-9 


van Zadelhoff et al. 2001 


JCMT 


01 


63//m 


2.7 X 10-*^ 


Thi et al. 2010 


Herschel PACS 


01 


145//m 


< 4 X 10"^ 


Thi et al. 2010 


Herschel PACS 


CII 


157/im 


< 4.8 X 10""^ 


Thi et al. 2010 


Herschel PACS 


Nell 


12.8 ^m 


3.9 - 4.7 X 10-^ 


Najita et al. 2010 


Spitzer IRS 






3.0 X 10"^ 


Pascucci &: Sterzik 2009 


VLT/VISIR (r < lOAU) 


Nelll 


15.5 iim 


2.0 X lO"'^ 


Najita et al. 2010 


Spitzer IRS 


H2 S(2) 


12.28 iim 


4.8 X lO"'^ 


Najita et al. 2010 


Spitzer IRS 






< 4.8 X lO"'^ 


Bitner et al. 2008 


Gemini/TEXES (r < 13.7AU) 


H2 S(l) 


17.0 fiui 


0.96 - 1.0 X 10^^ 


Najita et al. 2010 


Spitzer IRS 






< 5.5 X 10"^ 


Bitner et al. 2008 


Gemini/TEXES (r < 20 AU) 


OH 


MIR- SH 


6.2 X 10"*^ 


Najita et al. 2010 


Spitzer IRS 


COtiii, 


NIR 


1.0 X 10"^ 


Salyk et al. 2007 


Keck II/Echelle 




2.12/xm 


7.9 X 10"^ 


Bary et al. 2003 


KPNO/Phoenix 


H2 fluo. 


UV 


1.9 X 10"^ 


Herczeg et al. 2002 


HST/STIS 


01 


6300 A 


1.0 X 10-5 


Suzan Edwards 


Keck/Echelle 


01 


5577 A 


1.4 X 10"6 


Suzan Edwards 


Keck/Echelle 


SII 


4069 A 


3.0 X 10"^ 


Pascucci et al. 2011 


FEROS 
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Table 2: Dust model parameters 



Optically Thin Inner Disk 



Inner radius 


0.06 AU 


Outer radius 


3.5 AU 


Dust disk mass 


2.4 X lO^Sg 


Surface Density S 


oc 


Min. grain size 


0.9/um 


Max. grain size 


2.0//m 


Dust Composition 


Glassy Pyroxene 


Optically Thick Outer Disk 


Inner radius 


3.5 AU 


Outer radius 


200 AU 


Dust disk mass (with a < 1mm) 


4.8 X lO^^g 


Surface Density S 


Composite (Fig.l) 


Min. grain size {amin) 


O.Ol/im 


Max. grain size (amax) 


1 mm 


Dust Composition 


Silicates+graphite 



Table 3: Model calculations for lines (sum of inner and outer disk contributions) that constrain the 
inner disk mass 



Line 


Obs.Data 


SMH 


SML 






Contribution 






(M, = lO-^Mo) 


{Mi = 1O-5M0) 


{Mi = 10-6 


Mo) 


from r > 4AU 




(Lo) 


(Lo) 


(Lo) 


(Lo) 




(Lo) 


COvib (P Series) 


10-5 


2.1 X 10-5 


1.8 X 10-5 


1.6 X 10" 


-5 




OH (prompt) 


1.5 X 10"5 


7.8 X 10-6 


4 X 10-6 


2.8 X 10" 


-6 


2.8 X 10-6 


01 6300A 


10-5 


1.5 X 10-5 


8.4 X 10-6 


8.9 X 10" 


-6 


6.1 X 10-6 


( Thermal) 




(3.1 X 10-6) 


(2.3 X 10-6) 


(2.8 X 10" 




(1.1 X 10-6) 


( Prompt) 




(1.2 X 10-5) 


(6.1 X 10-6) 


(6.1 X 10" 




(5 X 10-6) 


01 5577A 


1.4 X 10-6 


1.92 X 10-6 


1.6 X 10-6 


9.4 X 10" 


-7 


7 X lO-'^ 


( Thermal) 




(2.2 X lO-'^) 


(1.0 X lO-'^) 


(1.3 X 10" 




(1.0 X 10-1°) 


(Prompt) 




(1.7 X 10-6) 


(1.5 X 10-6) 


(8.1 X 10" 




(7 X 10-^) 


H2 S(2) 


4.8 X lO-'^ 


6.8 X lO-'^ 


3.4 X lO-'^ 


2.5 X 10" 


-7 


2.2 X 10-"^ 


H2 S(l) 


1.0 X 10-6 


1.3 X 10-6 


1.2 X 10-6 


1.1 X 10" 


-6 


1.1 X 10-6 


Nell 


5.0 X 10-6 


3.1 X 10-6 


3.0 X 10-6 


3.4 X 10" 


-6 


2.3 X 10-6 


Nelll 


2.0 X lO-'^ 


3.0 X lO-'^ 


3.0 X 10-^ 


3.5 X 10" 


-7 


1.0 X 10-"^ 


01 linewidth (km s"^) 


10.0 


11.0 


6.0 


6.1 






a (vise, parameter) 




0.01-0.02 


0.1-0.2 


1-2 
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Table 4: Calculated model line luminosities for lines originating beyond 4AU 



Line (A) 


Obs. Data 


Model 




(Lo) 


(Lg) 


OI(63/zm) 


2.7 X 10"^ 


3 - 6 X 10-^ 


OI(145//m) 


< 4 X lO"'^ 


1 - 3 X lO-'^ 


CII(158//m) 


< 4.5 X 10~6 


3 X 10-6 


CO 6-5(430Aim) 


4.6 X lO-*^ 


8 X 10-^ 


CO 3-2(866/im) 


2.4 X 10"^ 


2.4 X 10-s 


CO 2-1(1. 3mm) 


7.9 X 10-9 


7.1 X lO-*' 


HCO+4 - 3(841/im) 


7.7 X 10"9 


4 X 10-9 
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Fig. 1. — The input spectrum used for TW Hya, which is a composite of the stellar, FUV and X-ray 
spectra shown here as the flux at the stellar surface as a function of photon energy. The vertical 
dashed lines indicate the EUV region which has been ignored in our modeling. 
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Fig. 2. — Surface density distribution of the SML (solid line) and SMH (dashed line) disks. The 
surface density follows a power law (S oc r~^) within 4 AU, rises over a "rim" region (3.5 — 4AU), 
and is approximately a power law (S oc r"'' '') from 4 — 40AU. In the outer disk (r > 40AU), the 
surface density rapidly decreases with radius. See §3.1 for the analytic formulation of S(r). 
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Fig. 3. — The density and temperature structure of the disk for the SML disk. The outer disk 
(f ^ 40AU) density structure shows a somewhat flattened disk structure, resulting in relatively low 
disk temperatures in these regions. The black line in both panels shows the height at which the 
visual extinction to the star, Ay = 1- The insets show the structure in the inner disk for r < lAU 
where most of the line emission is produced. 
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Fig. 4. — The vertical temperature(solid line) and density(dashed line) distribution is shown at 
a typical radius in the inner disk (O.lAU). The dominant heating (above) and cooling (below) 
mechanisms are shown and the radial extent where they are important in overall thermal balance 
is indicated. The photodissociation fronts of CO and H2 are also marked on the plots. 
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Fig. 5. — The disk structure in the mid-disk region at a radius of 6 AU. Here, the dominant FUV 
heating by PAHs and photoreactions near the Ay = 1 layer, and X-rays at the surface. Near the 
midplane, dust coUisions are important. The main gas coolants are CO and H2 rotational lines 
deeper in the disk, and 01 (63/im and 6300A) and [Nell] and [Aril] at the disk surface. 
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Fig. 6. — The disk structure at a representative radius in the outer disk (30AU). The dominant 
heating and coohng mechanisms and their radial extent is indicated in the figure. 
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r (AU) 

Fig. 7. — The cumulative line luminosities {L{r)) normalized to their total values (Table [6|) for the 
SML as a function of disk radius. The molecular lines, including prompt emission, usually arise 
from near the Ay = 1 layer, whereas the 01 and Nell thermal emission originate above this layer 
in the disk surface. The 10% and 90% levels (dotted lines) mark the radii between which most of 
the emission originates. For example, the [NeII]12.8/um line mainly arises from gas located between 
1 and 8 AU. 
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10 20 30 

J 

Fig. 8. — Calculated CO rovibrational fluxes for the v = 1 — 0, P(J= 1 — 30) transitions and the 
SML inner disk mass (solid line) are shown for each rotational quantum number J. The filled circles 
are the photosphere and veiling corrected measurements from Salyk et al. (2007). Open squares 
are corrected data from Rettig et al. (2004) (see text). The dashed line is for a higher disk mass 
(1O~^M0) that also fits the data reasonably well. The 1O~^M0 inner disk mass model (dotted line) 
is a poor fit to the data. 
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Fig. 9. — Contours of the thermal hne ratio OI6300A/OI5577A as a function of the gas temperature 
and electron density from a statistical equilibrium calculation of the O atomic levels are shown. TW 
Hya has an observed ratio of 7 (solid line) which requires high electron densities at typical neutral 
gas temperatures. For an electron abundance Xe ~ 0.01 and T ~ 8000K, gas densities n ~ 109 cm 
are needed to reproduce the observed 6300A/5577A line ratio by thermal 01 emission. 



